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Abstract 

A coherent state path integral of anti-commuting fields is considered for a two-band, semiconductor-related 
solid which is driven by a ultrashort, classical laser field. We describe the generation of exciton quasi-particles 
from the driving laser field as anomalous pairings of the fundamental, fermionic fields. This gives rise to 
Hubbard-Stratonovich transformations from the quartic, fermionic interaction to various Gaussian terms of 
self-energy matrices; the latter self-energy matrices are solely coupled to bilinear terms of anomalous-doubled, 
anti-commuting fields which are subsequently removed by integration and which create the determinant with 
the one-particle operator and the prevailing self-energy. We accomplish path integrals of even-valued self- 
energy matrices with Euclidean integration measure where three cases of increasing complexity are classified 
(scalar self-energy variable, density-related self-energy matrix and also a self-energy including anomalous dou- 
bled terms). According to the driving, anomalous-doubled Hamiltonian part, we also specify the case of a SSB 
with 'hinge' fields which factorizes the total self-energy matrix by a coset decomposition into density-related, 
block diagonal self-energy matrices of a background functional and into coset matrices with off-diagonal block 
generators for the anomalous pairings of fermions. In particular we investigate the transformation from the 
coset fields of a curved coset space, as the independent field degrees of freedom, to locally 'flat' fields with 
Euclidean integration measure. This allows to reduce the final path integral to solely 'Nambu'-doubled fields 
after a saddle point approximation for the density-related self-energy matrices and also allows to derive clas- 
sical field equations for exciton quasi-particles from various kinds of gradient expansions of the determinant. 
Despite classical evolution equations for 'Nambu'-doubled terms, one thus incorporates a semi-classical notion 
with quantum properties, due to the transformation and reduction to self-energy matrices which represent 
the irreducible parts of a diagrammatic propagation. The derived equations with various kinds of classical 
self-energy matrices also allow to examine self-induced transparency effects and the area theorem for ultra- 
short coherent transients of matter in combination with holography and selective Fourier optics. 
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1 Introduction 

1.1 Hubbard-Stratonovich transformation of quartic fields to self-energies 

Phenomena of many body physics can be mapped to different kinds of theories which can start from very direct 
approaches as in computational physics [I]-[l] or which can be based on various approximations. Apart from 
the YBE-related, integrable systems [5], there is a general lack of constant integrals of motion so that simple 
solutions of functions do usually not appear in many particle systems. In consequence, one has to rely on 
approximations, as by choosing a dependence of the self-energy on the irreducible diagrams or Green function 
parts or as by taking various manners of correlated field operator terms. The chosen approximation of the 
self-energy or correlated operators has a satisfactory, correct form if the prevailing physical phenomenon is 
considered within its computation; therefore, different manners of calculations can occur, each approximation 
being regarded as a particular theory aside from the original, beginning formulation with a fundamental Hamil- 
tonian, Green function and Dyson equation or the original definition of the correlations with the '/u// ' operators 
[6]- [15]. We prefer those "approximating theories" as the more appealing kind in many particle physics where 
one only reduces to a few terms, but still regards the essentials of the corresponding phenomenon. 
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In the present paper we investigate a coherent state path integral of fermions or anti-commuting fields, by 
way of example for a two-band, semiconductor-related solid which is driven by an external, classical electric 
field [6l [TBI HT] . We apply various kinds of Gaussian transformations from the quartic interaction of Fermi fields 
to various self-energy matrices. The latter operation or Hubbard-Stratonovich transformation ('HST') maps 
the quartic Fermi-fields of the original coherent state path integral to a bilinear term of anti-commuting fields 
with linear coupling to a self-energy matrix and a separate Gaussian factor of the corresponding self-energy. In 
section [2] we define the total Hamiltonian from second quantized, fermionic operators in normal ordering whose 
time development is determined by a coherent state path integral }16[ ITT] . Since the driving interaction with the 
classical electric field creates electron- hole pairs, one has to perform an anomalous doubling of anti-commuting 
fields for the presence of exciton quasi-particles within the original coherent state path integral. In section 
[3] subsequent HSTs lead to quadratic self-energy terms and to the linear coupling of the anomalous-doubled, 
bilinear Fermi fields with the self-energy. Three different kinds of HSTs are regarded within sections 13 . II to [331 
which we abbreviate symbolically in advance : 

Section [3TT] : Auxiliary variables o"(Tj''''\x) as scalar, real-valued self-energy (1-1) 



Section [3^ : Hermitian self-energy matrix S,, g/.^ JX,^^ ,X2;'J'j/^ , xi) of densities (1-2) 

Xu,A'^:?'\x2)(S)X;A'^-7'\x,) 'HST'secJionO S,,,.^,, (T^^^^^ \ f 2; T;:« ' , Xi ) (cf. ([HSiEISS])) ; 

Section [3T3l : Self-energy matrix S^^"^/.^ ..(t]^^^', afi) for densities and anomalous pairs (1-3) 



'HST' section [Ql ^ha (t^'^H'^ - . t'''ji ^ 



In section 13.11 we only transform to scalar self-energy variables with a contour time and spatial dependence, 
but without regard of any spin or electron-hole band labels. In section W% we extend the HST to a self-energy 
matrix which replaces dyadic products of Fermi fields for pure density-related terms of field operators. Section 
contains self-energy matrices, which follow from dyadic products of anomalous- doubled, anti-commuting 



fields and which thus consist of density-related, block diagonal matrix parts ^(T^^^^', ^2; 'Jj^^^'; afi) and 



also off-diagonal matrix parts ^^^"^ ^ ('^^2^^ ^X2]'J^_^^'^ ,xi); the latter substitute anomalous-doubled pairings 



of Fermi operators |16j-[19j. After the various HSTs in sections 13.11 to [331 and the removal of bilinear Fermi 
fields by integration, we obtain path integrals with the corresponding self-energy variables and matrices as 
the independent field degrees of freedom whose first order and second order variation straightforwardly result 
in saddle point equations and corresponding solutions for fiuctuation terms. We emphasize the precise time 
step order of the time evolution for a normal ordered Hamiltonian throughout the presented article and follow 
along the lines of the article |2Qj for a proper time step separation of field operators within coherent state path 
integrals. Section |4] has a more profound kind of HST where one transforms separately density-related and 
anomalous-doubled terms to corresponding self-energy matrices in combination with a spontaneous symmetry 
breaking ('SSB') and with 'hinge' fields [T6]; this leads to a factorization of the total self-energy matrix into 
density-related, block diagonal self-energy matrices and into anomalous-doubled, off-diagonal blocks. The 
latter are taken into account by coset matrices within a coset decomposition S0(9T, OT) / U(9T) ® U(OT) of the 
total self-energy matrix which itself has the detailed structure of a so(9T, 9T) generator within the orthogonal 
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group SO(ai,ai) (dimension 01 = (/i = 'e','h') x (s =t,|) x (r/j = ±) x (TV = t^/At) x {M^ = {L/^xf) = 
2x2x2xNx Afx', a factor of two for the electron- hole, spin and contour time metric degrees of freedom, 
respectively, times the number of time-like and discrete, spatial points defines the relevant dimension within 
S0(9T, 9T) /U(9T) ® U(*Jt) ). The coset decomposition combined with the factorization of the total self-energy 
allows for a projection onto the anomalous-doubled field degrees of freedom with the coset matrices whose path 
integral consists of the density-related self-energy as a background functional or as a solution from a saddle point 
equation. The coset decomposition involves various transformations, as the appropriate, invariant integration 
measure and as the transformation to Euclidean integration variables for first, second or higher order variations 
with classical coset matrices [21j. These involved appearing, but straightforward transformations are described 
in sections H] to [5] and are important for the appropriate consideration of the invariant integration measure, 
following from the square root of the coset metric tensor \16\ I21j . 

1.2 The importance of the appropriate, invariant integration measure 

In section H] we attain a coherent state path integral whose remaining field degrees of freedom, adapted by coset 
matrices for exciton-related parts, possess a nontrivial, non-Euclidean integration measure determined by the 
coset metric tensor of SO(0T, OT) /U(9T) (8) U(9T). In analogy we illustrate this problem for a multidimensional 
integral Z[x\ (ll.4atil.5cl) . (x = {x^, . . . , x^}), with an action A[x] where the Euclidean integration measure d[x] 
is modified by the square root of a metric tensor gij{x), similarly to the square root of the coset metric tensor 
in section a (cf. eqs. (I4.11ati4.12ll ) [22l[23] 

Z[x] = j d[x\ ^det[g{x)) exp{^yl[x]} ; (1.4a) 
d\x\ : Euclidean integration measure d[x] = dx^ A ... Adx^ ■ (1.4b) 



9{x) 



dx* gij{x) dx^ ; 



9ij[x) 

g'Kx) 



metric tensor with covariant indices ; 

inverted metric tensor with contravariant indices . 



(1.5a) 
(1.5b) 
(1.5c) 



The transformation to Euclidean variables dy^ = dyj (|1.6til.8p is related to the inverse square root of the metric 
tensor ^~^/^(x) as the appropriate Jacobi matrix where the symmetric property of the metric tensor allows 
a decomposition into orthogonal matrices Ojj(x) with det(0(x) ) = 1 and real eigenvalues A'^(x) {i,j,k,l = 
1,...,N) 



ids)'- 



dx' gij{x) dx^ = dx" Ofj,{x) A (x) Okj{x) dx^ = 

dx' (o^(x) • X^/^{x)^ (a1/2(x) • 0(^))^. dx^ = dy^ dy^ = dy^ dy 



(1.6) 



dy^ = dyj 
dji{x) dx" 



A^/2(x) •d(x))^ dx' ■ 
(\-^l\x)-dyy- ^y^=y\x) 



X 



x\y) . 



(1.7) 
(1.8) 



One obtains from (jl.Sp an arbitrariness of further orthogonal transformations, e.g. to dy'^ (|1.9ap . which 
are given by locally rotated Euclidean differentials, maintaining the local volume elements ()1.9bp . These 
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transformations correspond to gauge transformations and to the gauge invariance of the coset metric integration 
measure, determined in section S] for the anomalous-doubled field degrees of freedom 



dy'^ = dyr = {g^/\x))^^ dx' = (6^(x) • >}l\x) ■ d{x)) dx' = Oj,(x) {\^l\x) ■ 0(x) 



dx' 



= Oji^{x) dy^ = djf.{x) dyk ; {dy'^ = dy'j ; locally rotated Euclidean differentials); 
dljf] = d[y\ ; (invariance of volume elements under 'gauge' transformations) . 



(1.9a) 



(1.9b) 



The transformation to Euclidean differentials dy (or to equivalently rotated versions dy') involves the additional 
Jacobi matrix J*^ = (dx^/dy^) = ( 0"^(x) • A~^/^(x) )*^ whose determinant leads in combination with the original 
square root of the metric tensor to Euclidean integration variables. The simplified, Euclidean integration 
measure is accompanied by a transformation of field variables within the action A[x] to A'[y\ = A[x{y)] which 
alters the dependence of the field degrees of freedom within the classical action of the exponential 



■J k 

det(A) 
Z[x 



Qyk 



det 



d^(x) • \-^/\x) 



6^{x) . \~'/\x))\] = det[r^/'(x)] = (det[5 



[X 



-1/2 



d[y\ det g-^^^{x) J det{g{x)) [i A[x{y)]] = Z'[ 



d[y\ exp{i a' 



=1 



A'\ 



A\x{y)] ■ Z'[y\ = Z[x (y)] ; d[y\ = dy' A . . . A dy 



N 



(1.10) 

(1.11) 

(1.12) 
(1.13) 



In the following we compare the two cases T and '11' of variations for classical equations and their quadratic 
fiuctuation terms 'without' (case T) and under inclusion of the transformation to Euclidean integration vari- 
ables (case '11'). Case T disregards the nontrivial integration measure (det(^(if) ) )^/^ (|1.4am.5c|) in the first 
and second order variations so that one can simply expand the action A[x]. We thus attain from the first 
order variation the classical equations dA[x]/dx^ = ('in a transferred sense') whose solution (or 'solutions') 
determine the quadratic fiuctuation term (]1.14dp for the prevailing 'field' configuration x = xq 



Case 'I' : Z[x] = J d[x] ^det{g{x)) exp{i A[x]} ; 



neglect Ydet(^(x)) and expand A[x] for x 
dA[x 



xq following from : 



d'^Ai 



dx^ dx^ 



dx^ 
A[x] 

X=XQ 



; (i = l,...,iV) ^ f 



A[xq\ + (x* 



5x* 



xq ; (extremum); 

1 , . d'^A[x] 



X=XQ 



4) 



dx^ dx^ 



4) + 



x=xo 



quadratic fiuctuation term . 



(1.14a) 

(1.14b) 
(1.14c) 

(1.14d) 



In case '11' with the Euclidean integration variables we can similarly expand the action A'[y\ around a classical 
solution yo up to second order, but have to take into account the chain rule compared to case T. The first 
order variation (ll.lSbl) contains the multiplicative Jacobi matrix, compared to (|1.14b[) of case T, so that the 
classical solution yo directly follows from the transformation yo = y{xQ) according to a necessary 'rank A^'- 
property of the Jacobian for proper transformations (except the Jacobi determinant vanishes for particular 
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points which causes separate ranges of integration and transformation intervals). Nevertheless, this implies a 
different, second order fluctuation term (|1.15dp which is modified by a derivative and two multiplications of 
Jacobi matrices compared to case T 



Case 'II' : Z[x{y)] = Z'[ 



dA'[y\ 



dx^ dA[x 



d[y\ exp{^yl'[^/]} = / d[y\ exp{z ; 
dA[x 



y=yo 
because j\ 



dy^ dx^ 



dx^ 



0; (i=l,...,N) 



x = xo = x{yo) 



(1.15a) 
(1.15b) 



(^^) ^ ^ has to be of rank N for a valid transformation to Euclidean coordinates 

V dj/'- / y=yo 



A'[y\=A'[yo] + iy'-yh: 



dA'[y\ 



dy^ 



y=yo 



+ - (/ 
2! ^ 



Qyk Qyl 



(1.15c) 



A[xm=A[xiyo)] + iy''-yo: 



Qyk 



dA[x] 



4(/ 



Qyk Qyl 



dA[x] 



fluctuation term 



I 

' y—y'o 



dx'' 



y=yo 
+ 



dx^ 



y=yo 
+ 



xo=x{yo) 



Qyk Qyl 



xo=x(yo) 

dA\x] 



(9x* 


d^A[x] 


dx^ 




dy^ 


y=yo 


£o=x(yo) 


y=yo^ 



(y'-y(,) + ... ; 



dy'- 



y=yo 
dA[x] 



dx^ 



dx^ 

Qyk 



y=ya 



dx^ 



xo=x{yo) 

+ {J'k)y: 



d'^Ai 



dx^ dx^ 



y=ya 
d'^Alx] 



dx^ 

r 



^o=^(j/o) 



--yo dx"- dxi 



xo=x(yo) 



(1.15d) 



y=yo 



xo=x{yo) 



--yo 



In summary we can state that the classical solution, following from the first order variation, changes under 
the transformation to Euclidean integration variables because this induces a transformation of coordinates 
t/o = y(xo) with the square root of the metric tensor. The second (and also higher n-th) order fluctuations are 
altered by multiplications and by derivatives of the inverse square root of the metric tensor due to the chain 
rule for the transformation x = x{y) aside from 'gauge' invariant transformations as in (jl.9alll.9bl) 



Case T 
Case 'II 



2nd order fluctuations: 

x] 



d'^A\-' 



5x* dx^ 



X=Xo 

l/2/-„-r^ \i 



(y) 



dy^ 



dA[x] 



y=yo 



dx^ 



+ [9 



-1/2 



xo=x{y(i) 



d'^A\x\ 



^ dx'^ dx^ 



xo=x{yo) 



my, 



(1.16a) 
(1.16b) 



In analogy we perform in section transformations to Euclidean field degrees of freedom for the remaining 
anomalous-doubled self-energy parts within the coset matrices. In a similar manner one does attain changes for 
the classical field equations, according to the transformation of the field solution yo = y{xo) and to the 'gauge' 
symmetry (jl.9aHl.9b]) in a 'transferred sense', but has to consider the inverse square root of the corresponding 
coset metric tensor in order to conclude for the appropriate second and higher order fiuctuation terms. A 
further possibility is given by a modified classical action, according to 



(1.17a) 
(1.17b) 



Z[x] = J d[x] exp{zA[x]} , 
A[x] = tr[ln(g(f))] +A[i 



1.3 A glimpse in advance for the case with locally Euclidean coset fields 
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so that one allows for a further treatment of the nontrivial integration measure where one has to include the 
variations for the 'trace-log' term of the metric tensor. This straightforward procedure of variations necessitates 
various kinds of Christoffel-symbols of the metric tensor and also results in corrections, both in first and higher 
order variations, compared to cases T and '11'. This straightforward kind of variation can be found in random 
log gases of random matrix theory |24[ [25] . 



1.3 A glimpse in advance for the case with locaUy Euchdean coset fields 

In sections 13. II to [3^31 we specify HSTs which range from the transformation with scalar, real- valued self-energy 
variables to total self-energy matrices with anomalous terms in additional, off-diagonal blocks. However, we 
extend the latter case with anomalous-doubled parts and perform a factorization of the self-energy matrix into 
density-related self-energy matrices of a background functional and into coset matrices having generators of 
off-diagonal, 'Nambu' parts. This factorization is taken by a coset decomposition for exciton-related quasi- 
particles. Although the various steps of section S] and [5] are straightforward, we briefly list the resulting path 
integral of 'Euclidean' coset matrices for an overview in advance with regard to the familiar kind of a nonlinear 
Sigma model on a coset space S0(9T, m) / U(9T) (g) U(9T) | 



Z[3] « / d[a2(T;>;^T]:i')] 3ft-i(5);t(5);?^ 



xexp<; Ij^xJ2Y1 J2 

X j=0 ^=e,h 
-1 ^ 



2N+1 s=t,i 



X Tr 

a,b 



In ( i - [d^.^^^^zY"' + p a[t-i(5), t(2:)] (#C[^, (s 



(1.18) 



The independent field degrees of freedom are given by the matrices a('5';'i') whose derivative or variational 
increments 'a' are given in terms of quaternion- valued matrix elements AdlrL (T'''^i' ; T*''^l' ) for off-diagonal 
elements and, aside from an additional phase factor (see second line in (jLlQaj) ). by the field increments 
Aais^ (Tj V ; V ) = a( |as=±i('J'j^l')| exp{z5s=-|-i('TjV)} ) for the particular elements along the main quaternion- 
diagonal with element ('r2)^g^i. Apart from the contour time 'J^^^~^^ (j = 0, . . . , 2N + 1) and space vector x 
dependence, we include the semiconductor-related, two-band index fj, = 'e', 'h' and spin s =t, i, S = 2s = +1, —1 
degrees of freedom for the independent, locally Euclidean, 'Nambu' fields Aa('5';'i'); the already listed abbre- 
viations '5','!', . . . denote the total set of external and internal labels as the space-contour-time with additional 
band and spin indices, respectively 



a5('5';'i') = -Aa ('5';'i') = jz^fk^,, Aa^,(S^i^7';^J^J+ 

k=0 



+ (T2)At5Mi exp <^ t 



2i5s(^:f)i 

sinh(2|Ss(T^:')|) 



J 



(last term := tJ;;!; = ^'^^l^Ms := = , S = 2s) ) ; 



Si?(T;j;T;j) = |5s(^J)l exp05s(T%')} ; 



.('7,-)n 



(1.19a) 



(1.19b) 



^The dimension ^1 is given hy VI = {fi= 'e','h') x (s =t, I) x {r]j = ±) x {N = tiv/At) x (A^^ = {L/Ax^) = 2 x 2 x 2 x N x J^^; 
a factor of two for the electron-hole, spin and contour time metric degrees of freedom, respectively, times the number of time-like 
and discrete, spatial points defines the relevant dimension 91 within SO(?)'l, / U(9T) ® U(9t) . 
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, , S=2s {x} / r~(2)*.rr('7j).rr('?j)w rr(2)/rr(')j).q-('7j)>, 

4sS(T::::;;T;y)] = n n ( ^..^^..)^^-^..^^..) 

s=\,ij=l,...,2N^ * 
3 /S=2s,S'=2s' {^1,2} 

^ n n n (^p* • = s' = s) & (ji = & (xi = X2)) 

k=0 \ s,s'=t,i ji,2=l,.--,2Af 



(1.19c) 



2 z 



In comparison to [ISl [T7] , [201 [21] , [32] , we introduce a generalized gradient term {d^^^_gs^Z{'3';'2') ) ()1.20p within 
the determinant of (|1.18p which consists of the adjoint operator action of the sum of a one-particle operator 
and a saddle point solution (s)f,'l/.„ „(T. ^% :7. ) of the density-related background self-energy 



^ln(-^£)^('3';'2'))('5';'l'): 



(exp{[ln(-tf{[:^,(5)]),...]^} - i^Z^'^^'Cs'; 



ba 

'2') ('5';'!') 
~t, 



(1.20) 



The coset matrix Z^"" {'5' ■,'!') with locally Euclidean, 'Nambu' field increments Aa('5';'i'), Aa ('5';'i') in the off- 
diagonal blocks 6 / a is outlined in (|1.21a|) with the block diagonal, density-related terms aV'^"('5'; 'i') ()1.21bp 
which depend on the locally Euclidean fields of exciton-related quasi-particles 



A£-C5Vl')=f"^"^''^''^ A«('^V1') Y 
\ Aa ('5';'!') Ay^'^{'5';'i') I . 

3 

a9^1('5Vi') = -aV^^'^{'5';'i') = -i ^^(rfefa) x 



fc=0 



tanh 



X exp 



(-1)* 



tanh 



iss,(T;3;)hias,(T;3;)i 



^ ' ^ ibr 



sinh(2fe(T;;^;)i) st;; 

tanh('^--^^^^^^^':'^--<^^^')+tanh('^--^^^^^^ 



(k) /jt'iis' .y^Ji^^ 

SS^*!*^ i5,^5' 



X exp <i d7 / 



sinh(2|Ss,(T;:j;)i) st;:^" 



S5S1\ jf5,X5' jl,Xl' 



6^,^, 6s,s, tanh(|Ss(T]:£^)|) \?isi7^i')\ a5s(T,^J) ; 



•('7j)> 



(1.21a) 
(1.21b) 



(last term (7^^ := J^^^l = 7'^%Ms := = si , S = 2s)) ; (matrix)t = (matrix)*-^ ; 
(cf. gauge fixing in eqs. (|5.15a|5.15bp and appendix lA]) . 

The derived equations with various kinds of approximations can be applied for the investigation of self-induced 
transparency and the area theorem in the case of ultrashort coherent transients propagating in matter with 
reduced absorption [26]- [30]; this can be further combined with holography and Fourier optics. 
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2 Representation of second quantized Hamiltonians with coherent states 
2.1 The coherent state path integral for the time development 

In this section we define the total Hamiltonian Htot. (V'^j V^i t) of second quantized Fermi operators il^i_i,s{xi) (|2.1ap 
and transform to a coherent state path integral (|2.15|) which is used for the four different HST's in sections 
13.11 to 13.31 and HI We approximate the solid by a single pair of semiconductor-related electron-hole bands 
which are regarded by the indices /j., the opposite charges Qe, Qh (|2.1bp and further spin degrees of freedom 
s,s' =t)i (|2.1cp with projection labels a^, to be applied for the dipole moments. The corresponding anti- 
commutators of Fermi operators ip^^s{xi), '4'^^ s'i^"^) listed in ()2.1ap with their various definitions (j2.1bl2.1c|) 
for the single '(e)lectron', the single '(h)ole' and the two spin degrees of freedom. The delta function (5^"^^^ 
denotes a Kronecker delta Sg2,xi weighted by the total number J\fx = [-^Y °f space points in 'd' dimensions 
with system length 'L' and spatial grid intervals 'Ax'. The total number Mx = (^)'^ of space points can be 
applied as a parameter of expansion for the 'HST' transformed generating functions, provided that the density- 
related background fields do not vary considerably on the spatial grid with Mx points in comparison to the 
'Nambu'-related terms. This parameter Mx = (^)'^ of total space points is therefore similar to the parameter 
of matrix dimensions within random matrix theories. Spatial sums X^^. . . = Jx^di^)'^ ■ ■ ■ are normalized by 
the system volume 'L'^' to dimensionless integrals, whose overall dimensions are determined by the physical 
dimension of the prevailing integrand, apart from the primed, unweighted case ... of pure summation over 
numbered, spatial points. Double or multiple appearance of indices or spacetime points involve the standard 
convention of summations, except for indices or spacetime labels set in parentheses 

5 '4'ly{x2)}+ = V ^ss' ; 5 ^-^^,^'(^2)}+ = {^l,si^l) , tl^lyix2)}+ = ; (2.1a) 

H, U = '{e)lectron', '(h)ole' ; q^, : = -1, gr/j = +1 ; (2.1b) 

s, s' =t,i ; Of,, (3fj_ : ae = +1, an = ; /3e = 0, f3h = +1 ■ (2.1c) 

The total Hamiltonian Htot. (V'^; V*; ^) (|2.2ap is given in terms of a bilinear, kinetic density part Hf{''PK''P) 
()2.2bp with general dispersion t^^s's{x2) '^^^^^j which can be approximated by relations (j2.3a|2.3bp . and also 
contains an anomalous part He{iI^\iP', t) (I2.2cp . which is driven by an electric field, coupled to a dipole moment 
Pu,s';fi,s', hence, one only creates electron-hole pairs in correspondence with charge conversation. Moreover, 
the general dipole moment pu,s';fj.,s should be constrained to the two parallel spin directions s = s' ="[,-1 
in order to regard the spin one-property of the electromagnetic field with the two possible spin components 
2s = S = 2s' = S' = {+1,-1}. We therefore summarize this electron-hole generation by a rather general, 
time-dependent, driving dipole moment i>v,s';n,six2, xi;t) ^ 5s's 'I)^-^{x2, xi;t) (j2.3cp with the specific dipole 
inter-band moment Pu,s';fj.,s{x2,xi]t) ^ Ss'spu;fj.{x2, xi;t), a spatially dependent band gap Eg{x) and an electric 
field of central frequency uj and envelope EQ{x;t). According to various kinds of possible dipole moments, 
a rather general kind is used with the type ^'D^^s';fj.,s{x2, xi;t)^ and projection labels fS^, ctf^ for inter-band 
transitions which may be approximated by a constant dipole moment in the case of very ultrashort pulses [!3l]. 
Apart from the Coulomb interaction Hff{''P\'4') (|2.2dp of quartic Fermi operators and Coulomb potential 
V{x2.,xi)^ we also consider a random potential u{x,t) with a Gaussian distribution (12.3d|2.3ep for disorder 
and noise in the semiconductor-related solid. After an ensemble average according to (j2.3dl2.3ep . one thus 
also achieves an anti-hermitian, dissipative, quartic interaction of Fermi fields with the spacetime variance 
f{x2,t2]Xi,ti) and energy scale 'uq' 

Htot.(^^V';^) = i^F(^^^) + ^i^(^^^;^) + i^FF(Vi^^) + ^D(V^^^;^) ; (2.2a) 
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tJ.,S,S 



HF{4,li;)= ^i^l^,ix2) t^,s's{x2)5%% 



a 



^1,2 



(2.2b) 
(2.2c) 



(2.2d) 
(2.2e) 



M~\ (X2) 



rn,, ^ X2,xi I 



me, mh> ; 

M^^,^(x2) ■ generalized 'mass' metric tensor similar as in general relativity theory ; 
i>u,s';fi,s{x2,xi;t) ^ - 5s's(pu■,^^{x2,Xl■,t) • ^o(^(^2 + ^i);^)) cos||(/ia; - Eg{^{x2 + xi)) • t| 

i^i',s';M,s(^2,^i;t) ~ -Ss's[PuAx2,Xi;t) • ^o(^(^2 + Xi);t) 



exp^f^{fuj - Eg{^{x2 + xi)) -tj 



rotating wave approximation 



Pu-Ax2,xi;t) ^ 13^ ^^"2^ Ph;e(^i;i) ; 

7[u{x,t)\ =N'j> exp{-^J2xi/ ^(.^2,t2) f~^{x2,t2;xi,ti) u{xi,ti)} ; 



u{xi,ti) u{x2,t2) = UQ fixi,ti;X2,t2) 



(2.3a) 
(2.3b) 

(2.3c) 



(2.3d) 
(2.3e) 



We further extend the total Hamiltonian Htot^ip^ ,ip;t) (|2.2ap by a source term '^', defined in ()2.4|2.9a|2.10ap 
for generating bilinear observables of Fermi fields, and introduce the corresponding coherent state path integral 
(j2.4p on a non-equilibrium time contour which is taken into account by a 'plus' branch (j2.5ap for forward 
propagation (lower line in ()2.4|) ) and a 'minus' branch (j2.5bp (upper line in ()2.4p ) for backward propagation. 
In order to include the changing sign in the exponent of the time-step development operators, a contour metric 
rjj p.7p is defined in correspondence with a contour time (12. 6p which has equivalent values = 



2iV-j 



(|2.8ap . due to the separation of the sign with metric r]^ = —rjiN-j, (j / -^)- We accomplish a path integral 
by using the anti-commuting coherent states ()2.8bp with '(e)lectron'-, '(h)ole' and spin degrees of freedom and 
further distinction of the prevailing plus-, minus-branch so that the labeling XiM,s{'^j ^ ,x) with an additional 
sign '(^j)' has to be applied for distinguishing the Grassmann-valued coherent states between their plus-, minus- 
branch of propagation 



Z[3] 



exp{-l AtH^^^.^(Vit^^.i^)| . ... .exp{-jrAtR^^^.0,4,;h)] 



(2.4) 



2.2 Anomalous doubling of the one-particle part of the Hamiltonian 
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Contour time : '+' branch tff > . . . > tj+i > (t^ = T,) > . . . > t^ > to ; (2.5a) 

Contour time : '-' branch tg < . . . < tAr-(j+i) < (^iv-j = "^n+j) < • • • < iAr-i< (t]v); (2.5b) 

'^2N — to ) '^2iV-l = ^1 i • • • ) 'J'iV+2 = tff_2 ;'J']v+l = ijv_l = tjy j'^'iV-l = ^AT-l j • • • j '^1 = ^1 j '-''o = ^0 i (2-6) 

T/2]v = -1 ;??2jv-i = -1 ; • • • ;??jv+2 = -1 ; r/^r+i = -i ; 77^^ = +i ;t/]v-i = +i ; . . . ;r/i = +1 ;?7o = +1 ; (2.7) 



T — T 

•' j — •'2N—j 1 



= -??2jv-j ; (o<j<Af-i) ; r/iv = +1 ; 



(2.8a) 



Due to the over-completeness relation of coherent states and the normal ordering of the total Hamiltonian, we 
finally acquire the defining coherent state path integral ()2.9a|2.9bp for the time development 



2N 

m = (0| []exp{- 1 Atr/,- ti^^^ .^'^ ,^-7^+(^^^_^y2 

2N 2N 

X 



Z[3] 



(2.9a) 
(2.9b) 



x,ii,s ^ j=0 j=l 



2N 



i=i 

2.2 Anomalous doubling of the one-particle part of the Hamiltonian 

As one performs the additional ensemble average (j2.3dl2.3ep . one attains the ensemble averaged coherent state 

path integral Z[3] (|2.10aj) with the two-body interaction potential V(x2 5 'Ji2^^^ ^i, "Jji^^') (!2.10bj) . containing 

a dissipative noise part, and density pairs bC^^/^ ,xi) ()2.10cp for the quartic interaction of fields. An anti- 
hermitian epsilon term (e_|_ > 0) is added for analytic, convergent properties of Green functions following after 
the integration over the Grassmann fields 



2N 



^[a] = / n n ^^^^^^^H?^^^ 

x,fj,,s j=0 

2N 



(2.10a) 



^ ^ 2N 2N X 

xexp (^x;,.(^^\x) x,A'^'-'\S)-Y.xlA'j'-'\x) Xm,.(t]'T^\^)) 

x,^l,s ^ i=0 3=1 ^ 

. fi,u,s,s' 2N 

xexpJ- ^ "^Tjj At xt,s'i'^j''\^2) Su^i (^£^,s'4^2) - <5s's XM,s('^]-r'\^i)+ 

Xl,2 j=l 
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^1,2 il,2 = l 
Ml. 2, SI. 2 2Af+l 



Sl,2 Jl,2=0 

V(f2,T;^'^^';fi,T;;^^') = r/,, <5,,j, (<5^,^,^,J F(f2,^i) " § ^ ^^o / (x2 , T,, ; i'l , J ; (ji, = 1, • • • , 2iV) ; 

(2.10b) 



.{-iji-i) 



(2.10c) 



Note, that we have also specified the source term for generating observables in the last line of (|2.10a|) . However, 
aside from the density terms of Hf{iP\''P), Hpf{"iI>\"iP), and for H£){ip\r[^;t), one has the anomalous operator 
HEii^\Tp;t) with pairwise creation and annihilation of electrons and holes so that an anomalous doubling 
has to be taken for the anti-commuting fields. This is illustrated in relations (|2.11al|2.11gP by introducing 
the anomalous doubled fields ..^(t]^^^', xi) ()2.11ap . We emphasize the time step shift from 'ji — 1' to 'ji ' 



between the upper Xmi,si C^ji-i ^\ 



Xl 



1) and lower component x*fii,si{'^fi^\^i 



2) of the anomalous- 



doubled field g^CJj^^^ ,xi) (j2.11aj) . due to the normal ordering of the Hamilton operators. Furthermore, 
we separately define extensions (j2.11bl2.11cp without anomalous doubling for the end points 'ji = 0' and 
'ji = 2N + 1' in order to consider the exact, proper sequence of time steps for the quantum problem (cf. Ref. 
|20]). This extension has also to be regarded for the hermitian conjugation in ()2.11dp with the peculiar 
time step between upper '6 = 1' and lower component '6 = 2', but in reversed order compared to ()2.11ap . The 
particular extension for the end points 'j2 = 0' and 'j2 = 2N + 1' is similarly given in Eqs. (|2.11el2.11fp . also 

with the absence of anomalous doubling as in (j2.11b|2.11cp . The hermitian conjugated field Hj'i'2,s2('Jj2^^\ ^2) 
()2.11dp also follows by multiplying S^^ ^^^t]^^^', xi) (|2.11ap with the 'Nambu'-related Pauli matrix (fi)''" and 
a subsequent transposition so that one has property (2.11g) which is used in removing the 'Nambu'-doubled 
anti-commuting fields by integration 



^a=l /'<j(''(ji=0)=- 



^/^l.siV {ii=0) 

a=l /q'(''(jl=2JV + l) = 

0i=2iV+l) 



""M2,S 



-tt,b=l/rr(''j2=0=- 
>2,S2\^''02=0) 



,Xl 
,Xi 

,X2 
,X2 



"M2,S2l •^(j2=2Ar+l) '-^2 



( XMi,si(^ii-i '\^l) 

V x;:,.,(^^\fi) 

^f^^M {ji=2N) '-^1 



( x;;2,«2(7'-2''',^2) , XA.2,s2(^.?-r'\^2) ) ; 

^fi2,S^-' j2=0 '• 



X2 

, ^(jlj2=2N = -) ^ 
Xn2,s\-' j2=2N '-^2 



(2.11a) 

(2.11b) 
(2.11c) 
(2.11d) 
(2.11e) 
(2.11f) 
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.(I,, ) 



(2.11g) 



As one applies the anti-commuting fields ..^(T)^'^^', xi) (j2.11ap . Hj'j'j.saC'^ja^^^ ^2) ^l^-Hdp to the exponent in 

(fZTOill . one obtains relation ([232]) with doubled one-part icle part H|^" ^,.^^^ (f 2 , 'J'j'a'' ' ; ^1 , T-j'^^ ' ) (|2.14all2.14d[) 
which has to be combined with the 'Nambu' metric tensors (I2.13al2.13b|i ^ for the process of anomalous doubling 
of anti-commuting fields (cf. Ref . [20j ) 



. 2N 2N 

+ Pua,i Xu,s'{7'^-i^\x2) I),,,s';/i,4^2,^i;'J'j+(,,^.-l)/2) X^.s C^j-l"'', ^l) + 

I ^/i2,S2V'^J2 1-^2) Ilj2 2 <^ t^2,S2■,^^l,s\\•' 32 )'^2)-^ji j-^lj '/ji ^/iijSil-Jji j-^l; 
^1,2 Jl,2=0 



(2.12) 



1 



tJ;iy,s,s' 2N+1 



^1,2 ii,2=o 



+ ^.2fc^M,s(^.2" >^2;T::-\xi)ry,, 



= diag{^ , ^} ; (2.13a) 

a=l a=2 

S'^' = dab diag{^ , ^} = • F ; (2.13b) 

a=l a=2 

+ ^ X '^i2ji % '5v/.(£^,s's(^2) - I Vj2 ^s's ' 

H^s';,,si^2,7^''^^^^'^'') = {f^s'^,,s{^2,7^''^^^^'^'')f = ^s's Ci ('^....1+1 - ^J^^n)+ (2.14b) 

+ * X ^j2,h ^v^, ( (£a*,s's(x2) - ^ e+) ! 

Hi^s';M,s(^2,T,f '^^1,T]^''') = ^ dj2,h Vj2 (^u ^Iy;/.,s(^2,Xi;Tj2 + (^^,^_l)/2) ; (2.14c) 

H^)^5';;.,.(^2,^,?'';^1,t]'''') = ^ '^i2,ii /^J^a^ i^!^,s';M,s(^2,Xi;^ji + (^,^-l)/2) • (2.14d) 

The combination of (j2.10ap with (j2. 12112. 14dp results into the path integral (j2.15p which contains the one- 
particle part in an anomalous-doubled manner and the remaining interaction part of quadratic, density related 
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pairs ()2.10cp . each composed of fermionic field pairs for density terms 



zm = y n n - %r^'°^ - ^ (2.15) 

x,/i,s i=o 

xexp<-- 2^ 2^ =i.,s'('^j2 '^2j [1 ,.(X2, ,2:1, ji + 



^1,2 il,2=0 



^ 2N 

Xl,2 il,2 = l 

3 The various kinds of Hubbard- Stratonovich transformations 

3.1 Auxiliary variables cr('jj^^\x) as scalar, real-valued self-energy 

The simplest possibility of a HST is given in (j3.1ap for transforming from densities b{7j^^'^ , ^2) • • • 1 ^1) 

or quartic interaction of anti-commuting fields to a bilinear term of fields with a single factor b{7j_^^'^ ,xi) 

and auxiliary, real variable ai7j^^^ ,xi). This is achieved by the additional Gaussian term (j3.1ap of the real 
self-energy variables whose 'variance' is determined in relation (|3.1b|l with the 'inverted', ensemble-averaged 
potential "^~'^{x2-,7^^^^^ ; xi, t]^^^'). The bilinear term of anti-commuting fields b{7^^^'^\x2) with self-energy 



variable a{7j^^^ ,xi) is converted to the anomalous-doubled form (jS.lcp with Grassmann variables ()2.11al 2.11g) 



and the 'Nambu' metric tensor S^"^ ()2.13bp 

2N 

Xl,2 il,2 = l 

= jd[<y{7f^\x)] exp{^^ E ^(^2'''' ^2) (x2, tJ^^^^^ fi, T^:^^') ^^^^ 

jl,2=l XI, 2 

2N 

xexp{-zAi ^ ^b(T];'-',f2)- 2 5,,,, r^,, • ^(7^:^', f 1)} ; 

ii,2=i 2 

2N 

^^V(x2,T;:-';f4,T;:-^) r?,, V-^{x,,7^^,x,,T^^') = 5,,,, r/,, 5f^% ; (3.1b) 
i4=i X4 

27V 

E E KtS-\x2) • 2 r?,, . ^(T]:-\fi) = (3.1c) 

Xl,2 il,2 = l 

^,v,s,s' 2N 

= E E ^2) 2 5,^5,,,4^^^5,,,,,r?,^ a(T;:^'^\xi) Xm,s(<V'^ 

Xl,2 il,2 = l 



3.1 Auxiliary variables (T(Tj''^',x) as scalar, real-valued self-energy 
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^1,2 il,2 = l 



As one applies the exponential form p.2p of ()3.1cp . one finally acquires the coherent state path integral p.3p 
with only bilinear, anomalous-doubled anti-commuting fields (|2.11al 2.11g ) instead of the quartic interaction of 
Grassmann fields as in (j2.15p 



(3.2) 



. 2N 

Xl,2 il,2 = l 

. . 2N . 

il,2 = l £'l,2 

^1,2 il,2 = l 



Z[a] = / d[a(T^',x)] exp<!^ 



2W 



^ E E^(^''^^2)r/,2v-^(f2,T^^^l,T-^o%.^(T^^ (3-3) 

jl,2 = l Xi_2 

J, ^^,v,s,s' 2N+1 

77; — exp{-l E E 

x,^,sj=0 a?i,2 ii,2=0 



27V 



X 



+ <5.;. S,,, rj,, 2 (zA*) 4^^^ 6,,,,, q, S"'^ 



-| ba 



The bilinear Grassmann term in (13. 3p is removed by the property of anti-commuting integrations where one 
has to use the identity p.4p for an anti-symmetric matrix M of even-integer dimensions as a one-particle part. 
One can also state that the symmetric part of the matrix M cancels in the exponent of (13. 4p so that one is 
only left with the anti-symmetric part of M in the resulting square root of the determinant (det[M])^/^ 



~ 2n 

J d[^,] exp{-i e,} = (det[M 



1/2 



M = -M^ . 



(3.4) 



As we define the matrix Mj.^ ^{7^^^^ ,X2;7j_^^^ ,xi) (|3.5ap and use property p.llgP between the hermitian 

conjugated fields Ef}l^s2{'^j2'^\^2) and their transposed form ( (fi)^'* H^2,s2('^j'2^^^' ^2) )^, we can reduce the path 
integral ([331) to ([3l5b]) with the square root of the determinant as exemplified in ()3.4p . Note that the Pauli 
matrix (fi)^" appears as a redundant factor within the determinant and can be disregarded due to the overall 



even dimensions of the matrix Afj^^^, (T)2^^',X2;T)j^^\xi) so that the path integral (l3.5bp is accomplished 



JVn-, ) 



with the remaining, real- and even- valued self-energy variables a{7j_^^^ ,xi) 
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(i[(T(Tj''^\ x)] exp 



+ <5.M Ss's Vn ^) Ci ^n^n ^(^.^ , ^i) ; (3.5a) 

r 1 

E^(<''^^2)r?,,V-i(f2,T;^^^^';xi,T]:^^'),7,^^ (3.5b) 

il,2 = l XI 2 



f DET 



1 A>(3) ^^^("^2) - .n-t-'ji) 



1/2 



^ a{7^T'\x) = x) . (3.5c) 

The first order variation of tlie patli integral ()3.5bp witli respect to (5cr(Tj^^^ ,xi) results into the saddle point 

equation (13.6P where one has to scale the self-energy variable fTo('T|^''^\ xi) by the value Nx ^''^ in order to 
achieve proper solutions according to the total number of spatial points 

2JV ^l,s 



i-(''j=2]V+l) 



CTo(T,^ ^xi) _ ^ 



Ml" 



J2 = li £2;(<'=1.2) 



An 



-l;aa 



7^'\x2^7'^'\x2) Vn S"'^ V(x2,T;^'^-^^xi,T;:^-^^) . 



(3.6) 



Solutions of (j3.6p follow from iteration with (To(X,/^ ,xi)/J\fx or continued fraction which involves the eigen- 
value problem with the operator 



2N-\-l tJ.,s 

EE 

ii=o XI 



(3.7) 



+ 6., 6,,, r?,, 2(z A*) 4^^^ 5^.^^^.^ s"^ ao(T]:«\ xi) 



= E'Af ^M u s'i'-^j2^^ 5^2) ; (eigenstate label M with complex eigenvalue .Em), 
so that the first order variational equation (j3.6p takes the form 

2N ^J.,s a=l,2 ,y,a {^^'-^32^ - 

4EE E 



J2 = l X2 Mgeigenstates 



{Em/Mx 



l/2^ 



(3.8) 



In general one could have problems to find proper solutions of (I3.7|3.8p : however, one can use a rotational or 
translational invariance in most cases so that a 2D or 3D problem simplifies to radial equations with angular 
momenta labels or also wave- vectors. The iterated terms of p.6|3.7p yield a complex solution aQ{7^^^\x) whose 
imaginary part has to comply with the sign of the imaginary (e > 0)-terms for a proper, stable propagation or 
the correct analytic convergence of Green functions. This implies a complex-valued eigenvalue problem (j3.7p 
whose imaginary eigenvalues Em have also to correspond to the sign of the infinitesimal, imaginary (e > 0)- 
term. If the imaginary parts of the solution (To('J, ^ , x) or Em take opposite sign as the imaginary {e > 0)-term 
of the Green function, one acquires an instable range of parameters for the semiconductor-related solid which 
may have a physical origin or may be related to insufficient approximations with our simplified HST and 
auxiliary, scalar self-energy variable cj(T,- ■' ,x). This range of instable properties (also due to possible artifacts 
of our HST in this section) can be improved by the more detailed kinds of HST's described in sections 13.21 to 



(v ■ ) (') ■ ) 
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3.2 Auxiliary matrix S^,,/;^ ^('J,^^^', 0:2; "^ji^^ , a^i) as hermitian self-energy density 

We extend the real self-energy variable a{7j \x), replacing the sum ^^^sXl^^sC^j ^ -.x) X^J.,s 

a self-energy matrix of fermionic densities which is specified as a dyadic product (|3.9ap of Xtii,si{'^^j[^-i^\xi) 

with the complex conjugate X^2,«2('^i2^^'' ^2)- The self-energy density ^y^s'■,^l,s{'^%^^\x2',7'^^-^^^\xl) ()3.9bp has a 
hermitian symmetry where the hermitian conjugation involves the electron- hole and spin labels with inclusion 
of the time contour and coordinate variables. As we apply the operation of dyadic products to the quartic 
interaction of Fermi fields, one attains the interaction ()3.10p in terms of two, solely density-related matrices 
from the dyadic product (|3.9aP , which allow for a similar HST as in section 13.11 but with matrices for densities 
in place of scalar variables representing sums of fermionic densities 

-Rmi,^i;/^2,s2(^]i''\^i;^.?'\^2) =Xm,si(^-V'''^i) ^ X%,s2{'^t'^^ ,X2) ; (3.9a) 
^/ii,si;/i2,s2 ('^ji ^ '^i!'^j2 ^ 5^2) = ( ^/ii,si;At2,s2 ('^ji ^ iXi;7j^^ ',X2j | ; (3.9b) 



2N 



{')„) 



Ati,2,si,2 2N 



Xl,2 il,2 = l 



^1,2 il,2 = l 



(3.10) 



^ Xm2,^2(^^2''''^2) % 9^.2 X/.2,S2(^]2-1 '''^2) X^i ,Si (^,'1'' ^l) ^/jl^m ,^1 (^il " 1 ' ' ' ^l) 

M,2,si,2 2N 



E E v(x2,^:-^;x„t;:-^ 



^1,2 ii,2=i 

One has to start from a similar Gaussian identity with a hermitian self-energy matrix of densities as in the case 
of scalar self-energy variables (t{7j ^ ,x) in section [3?T| we then transform the quartic interaction of Grassmann 
fields to bilinear fields given as a hermitian dyadic product with linear coupling to the hermitian self-energy 
matrix in a trace relation of internal and time contour, coordinate variables 



exp{-^^E E K^^<-2)^2^..K<"'>-i)^(-2,T;:-';fi,T;:-')} 

^l,2il,2=l 



(3.11) 



Ml, 2, SI, 2 27V 

d[%,s'-,^A'^l''\x2-,7f:'\s^)^ xexv\-if E E ^~'(^2,Tj^';fi,Tj^') X 



Sl,2 il,2 = l 



^32 ^A'2,S2;A»1,S1 ( '^j2 ^ ^X2','7j^ ^ jXlj Tjj^ Qf^-^ S^j ,5j;^2,S2 ( '^Jl ^ :Xi;7j^^ ,X2 



Ml, 2, SI, 2 2N 

xexp^2f E E 

^1,2 jl,2 = l 



'/j2 ^^2 ^ V2,S2;mi,si I -J j2 ' 2i •Jjj^ j-^ll 'hi y^i ^Mi,si;/i2,S2 I ji '■^1''^J2 ' 2 



Since the one-particle part H|^'*s/.„ s{x2,7j'^^^'' ; xi,'j]^''^') ()2.14all2.14dp also contains anomalous terms with dipole 



moments, one has to perform the anomalous doubling within the coupling of " (dyadic product) x (hermitian self- 

"'!^,s';/j,sV j2 



energy matrix)"; hence, we have to introduce two hermitian self-energy matrices 5]^"^, ('J^/'^ ,X2',7j ''^ ,xi), 
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3 THE VARIOUS KINDS OF HUBBARD-STRATONOVICH TRANSFORMATIONS 



(a = 1, 2) which are related by transposition with opposite sign (j3.12a|3.12bp . One thus achieves relation (|3.12cp 



with the anomalous-doubled Fermi fields ^^'^/(t]^''^', X2) 



(v ■ ) 

('J'j/^ ,3^1) and the two self-energy matrices 



which only comprise density terms without any 'Nambu' part 



t 



f.22 /rrt-'ia) ^^.V^n) - ^ 



exp^2f 



E E 

a;i,2 ii,2=i 



??j2 9/^2 -^M2,s2;Aii,si ( '^^2 '•^2;'J'ji ) j X 



(3.12a) 
(3.12b) 
(3.12c) 



^ ^Ji 9aii ^Mi.si;/i2,S2 ( "^Ji ^ ^ ,X2 



Ml, 2, Si, 2 2Af 

exp^-.2^ E E X 

^1,2 il,2 = l 

^ fj.,u,s,s' 2N 

X 9m S^i,si;/.2,s2('J'ji''\^i;'J'J'\^2) ??,2 Qm2 XM2j2('J'J-r'^^2) ^ = expj -i ^ ^ 

^1,2 il,2 = l 

X Htt';,,(T];'-\x2) ry,, g^, 2 (. f ) x^; T;:-\ fi) r?,^ H^,,(t;:-\ fi) } . 



After insertion of p.lip . (l3.12cp into ()2.15p . we obtain the path integral ()3.13aP whose 'Nambu'-doubled, anti- 
commuting fields are removed by integration according to the already given relation (j3.4p in section 13.11 The 
resulting path integral ()3.13bp only contains the hermitian self-energy density Sj^^^/.^ ,,('J'j2^^\ 2^2! 'Jj^^^^ ^i) as 
remaining quantum mechanical field degrees of freedom 



Z[d] = / d[E,,,.^,,(T];'^-^',f2;^-:^^',xi)] exp 



Ml, 2, SI, 2 2N 



At 



E E V-i(f2,<^^';xi,T];'^-^^) X (3.13a) 
^■1,2 ii,2=i 



'/j2 ^^2 ^At2,s2;Mi,si \ j2 i^'i'-'n '/ii y^i ^Mi,si;M2,s2 I ji j-^lj-Jja '-^2 



2Ar 



lb /rr(''j2^ 
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J2 



X2 X 



i'*H^';;.,s(^2,T;-^';xi,T;;-')I- + r^, 0;.,.(^2"^^2;T;;-^Xl) r?,, + r?,, 2^^f 

-1 fea 



i6a 



(t],^-\x2;^^'^',^i) H^,,(T;:-\xi)^; 



Z[3]= / d[S,y;^,,(T2^^^',f2;T];^'^\fi)] exp 



Mi,2,si,2 2N 



f E E V-^f2,T;:-';x„T;- , X 



(3.13b) 



^J2 ^M2 ^M2,S2;M1,S1 ( '^32^^ '^2;'J'j/^ , 2^1 ) (Z/Xi S 



^1,2 il,2 = l 
^M1>Si;M2,S2 I -^31 '•^l'"'j2 ''''2 
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ha \ 1/2 

- u,s';iJ,,s J 
■(jl'or'j2 = 27V + l) ,(cf. (in 

The final path integral ()3.13bp consists of the matrix -^//^l!'^^/.^ s('Jj2^^\ 2^2; "Jj^^^^ a^i) P-14p whose anomalous 



X2 ; \ f 1 ) EE for ( ji 'or'j2 = 0) 'or' ( ji 'or'Ja = 27V + 1) , (cf. (12. 1 lbl2. 1 If p ) . (3. 13c) 



terms (6 ^ a) are still restricted to those of the one-particle part H^^"^ ^(^2) 'J'jV^^'; In sections 13.31 



[4] we generalize to self-energy matrices which also comprise the anomalous terms ^(T*^^^', X2; 'Jji'^^ , 

in the off-diagonal 'Nambu' blocks (6 ^ a) 

Lj6a 



(3.14) 



In analogy to section ISTTl we take the first order variation of (|3.13bp so that one achieves the saddle point equation 



(I3.15ap in terms of the self-energy density S^^]"^ ^(T]J''^\ X2; 'Jji^^ 'j 2^1) whose solution is also determined by 
continued fraction and a similar eigenvalue problem (j3.15b|3.15cp of complex-valued matrices in compliance 
with the imaginary [e > 0)-terms of stable propagating Green functions (relations p.l5a|3.15cp are free of any 
summations over spacetime coordinates and internal state labels, aside from relation ()3.15bp with the weighted 
spatial sum . . . = Yl'xi ■ ■ ■ ^^"^ contour time step summations) 



V-i(f2,^-';fi,T;:-') Si^S....(^"^^2;<-\i:0/A/; 
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//;/i2,S2;Ml,Sl ^ J; 
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ji ,Xl) ^^^II-fj,j^,si;ti2,S2^ h '^1'''J2 '^2j 



2Af+l M.s 

ii=o 

(eigenstate label M with complex eigenvalue i?Af) ; 



Em Ku,M-7'\S2 



M2,S2;M1 



MSeigenstates 



(a=l) 



(a=2) 



(3.15a) 



(3.15b) 



(3.15c) 



-T,(a=2).^('7j2) 



M;/i2,S2('^J2 ' 



,X2] 



Em 



Em 



3.3 Self-energy ^(Tj^^^', X2; T^^^^ , xi) of densities '6 = a' and of anomalous parts '6 7^ a' 



Since the one-particle part H^°g,.^ ^(5:2) "J^J^^ ; 2^1, 'J'j/^ ) contains anomalous terms, one has also to introduce 
'Nambu' parts in the off-diagonal bloks (6 ^ a) of self-energy matrices. We therefore take the dyadic product 

(I3.16j) of the anomalous-doubled, anti-commuting fields 

anomalous-doubled, fermionic density pairs 5S(T^V''' ,xi) (ISTTP with extended 'Nambu' metric tensor (f3lB 
which implies an overall, additional factor ^. The combination of (|3.16H3.18]) transforms the quartic interaction 
of Fermi fields to a trace relation (j3.19p of density pairs 3?^^^ .,^.^2,s2('^j*i^^'' ^i' '^j'2^^'' ^2) (|3.16p with inclusion of 



'-'1ii,si{'^fi^\^i) and 'E^ill^s2{7^^2'^\x2) and consider 

ba 
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3 THE VARIOUS KINDS OF HUBBARD-STRATONOVICH TRANSFORMATIONS 



'Nambu' parts which involve the additional trace summation (a, b = 1,2) 



J^n ) 



Xl,fl,S 



aba aba 



2N 



2N 
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Xl,2 Jl,2 = l 
Ml, 2,^*1, 2 2Af 



^1,2 il,2 = l 



X Tr 
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^M2 'J2 ■^/X2,S2;M1>«1 J2 )-^2,-^ji ' ■^1'' "^Ml 'Ol ^M1>Si;M2,S2V Jl '■^ll'^J2 '"^S^ 



(3.16) 



(3.17) 
(3.18) 



y: e ,,2 vn K^"^^l) v(f2,<-';fi,Tj:-^) = 5; vix.^^^^x.^'jf^') x (3.19) 



According to the dyadic product (13.160 of anomalous-doubled, anti-commuting fields, we take into account 
'Nambu' parts (6 ^ a) of an overall, hermitian self-energy matrix SJJ^ ,,^.^2,s2('^ii^^^'^i''^J2^^''^2) (|3.20ap whose 
diagonal blocks (13.20bl3.20cP are related by opposite sign and transposition and whose two, anti-symmetric 
'Nambu' parts are related by hermitian conjugation (j3.20dl3. 20^3.201]) 



yab (cj(^h'> ^_,.t'''^2) 

^/ii,si;/i2,S2V ii '■^l''^J2 ' 1 



/ii,si;/i2,S2 V ii 
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J1„ ) 



/il,si;/i2,S2 V Jl 



J2 



yll /n-f^ii) ^^.rr(''i2) ^\ yl2 /n-'^ii) .'T*''J2) ^^ 



"/il,si;/i2,S2 V Ji 
^21 



Mi,«i;/i2,s2 ' _ 

T-^^n) - .<T(''i2) - \ v22 /o-f^Ji) .rr(''j2' - \ 



/J1iSi;M2,S2 



(t^:-\x,;t;:-\x2))'; 



^A»i,si;At2,S2V ii '■^i''^J2 ' 
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yl2 .^"^2) 



y21 /rj-f^n) ^,.T-*''^2' 



(3.20a) 
(3.20b) 

(3.20c) 

(3.20d) 

(3.20e) 

(3.20f) 



(for (ii'or'j2 = 0) 'or' (ji'or'j2 = 2iV + 1) ) , (cf. (|2. 1 lbl2. 1 If M . (3.20g) 



As one applies the dyadic product ^5ii,si;/i2,s2('^ii^^'' ^i' '^j'2^^'' ^2) (|3. 16113. 15|) and the self-energy matrix (|3.20at - 
I3.20fp with the 'Nambu' parts in the off-diagonal blocks (a ^ b) to the original, quartic interaction of Fermi 
fields, one attains the HST (|3.2ip with Gaussian term of self-energy matrices and a phase factor with linear 
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coupling between dyadic product of doubled Fermi fields and the anomalous-doubled self-energy in the exponent 

2N 

. '.'y '/y . ^ 'y t ^ \ (3-21) 

^1,2 Jl,2 = 
Ml, 2, SI, 2 2N 

I , A/ 
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^ ^Ul %1 
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M1i'5i;M2,S2 \ " Jl 



■^^2 ^J2 

Mi,2,si,2 2Af 

xexpiff E E 

^1,2 il,2 = l 



M2,S2;M1,S1 \ J2 '■^^i-'ji ^ 



X Tr 
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^U2 %2 ^ 



ba 



Ai2,S2;Mii'5i V ^^2 



Sa6 

Vi,Si;M2,S2 V " Jl 



('',1 ) 



.2 '^2 



Insertion of the HST p.2ip into the path integral ()2.15p results into relation p.22ap with additional Gaussian 
factor of anomalous-doubled self-energy matrices, but with only bilinear, anti-commuting fields. These are 
removed by integration according to relation p.4|) so that we finally achieve path integral (I3.22bp with the 
square root of the determinant and the 'Nambu'-doubled self-energy matrix as remaining field degree of freedom 



Z[3] = /4E^>^^,,(T*^^-',X2;T;^-\xi)] exp 



Ml, 2, Si, 2 2N 



E E V-^(x2,<-^x„T;:-') X (3.22a) 
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X exp 
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\ E E -K'(^:-\^2) 

^^1,2 il,2=0 



-1 ba 



z[a]= /4i:^>^,X'" '^2;t]:«',xi)] expj-1^ E E ^"'(^2,t;:-^xi,t]: 

^1,2 Jl,2 = l 



Ml, 2, Si, 2 2N 
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4 COSET DECOMPOSITION WITH DENSITY- AND EXCITON-RELATED PARTS 



The saddle point equation of (|3.22bp and its related eigenvalue problem follow in analogy to sections I3.H 13.21 
and are therefore listed in brevity in subsequent equations (|3.23ati3.23fp (relations (j3.23b|3.23c|3.23e|3.23f|) 
are free of any summations over space-contour-time coordinates and internal state labels; note, however, that 
()3.23dp implies the scaled spatial and contour time step summations) 



T{V(a)6a ^^^("^2' - .q-^^n^ - N _ jbb Ciba ( ^("^2) - n-'^Jih faa , „ ^ba /rr(''j2' - . n-t^Ji > -f ^ 
_L Qbb At\ yba To-T'"-'!' T 1 S 

V-'^l^ n-'^M* - rr(''jl)x v(0)ll /q-f-'M' - . 
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t];'"',xi)/aa,= 
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(eigenstate label M with complex eigenvalue Em) ', 

^ ^^:Lc^^'\-2) ^Z:^cj^^\x^ 
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Em 
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E 

M Seigenstates 



;A'2,^2'- J2 '-^2; ^M;fii,si\-'3i ^Xl) ^ ^M;/ii,siy-'ji '■^l/' ^ M;H2,S2^ ^' 

Em 
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^ M;H2,S2^'' ^2 '-^2; 



mi + 

(3.23a) 
(3.23b) 

(3.23c) 

(3.23d) 
(3.23e) 

1 

(3.23f) 



4 Coset decomposition with density- and exciton-related parts 

4.1 SSB with 'hinge' fields and anomalous-doubled parts 

In section [3T3I the density and 'Nambu' parts are regarded within a single, total matrix, having the appropriate 
block structure (I3.20ap . In this section we aim at a further factorization of the total self-energy into block 
diagonal, density-related parts and coset matrices whose generators are composed of the anomalous-doubled 
field degrees of freedom. In order to separate the coset matrices with 'Nambu' generators from the densities in 
a HST, one has to introduce block diagonal self-energy densities s°2,s2;Ati,si ('^j'2^^'' ^2; ifi) (I4.1a|4.1bp in a 

SSB as invariant ground or vacuum state and a total self-energy matrix .j^.^^ .,1 ('^J2^^'' ^2; 2^1) (|4.1cp 
with anti-hermitian 'Nambu' parts for the factorization with coset matrices (cf. section 3 in [16j). The ma- 

trix <^S^"2.s2;Mi.^i(^t' '^^2; T;;^"^\xi) dM has a similar structure glHlFTiD as E^'^^,,^.^^,,^(T;^'^^\x2;T];^^\xi) 
(|3.20ap in relations (]3.20bll3.20fp of section 13.31 except from the inclusion of the imaginary factor 'z' in the off- 
diagonal blocks for the anti-hermitian property which is required for a factorization by a coset decomposition 



4.1 SSB with 'hinge' fields and anomalous-doubled parts 
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(4.1b) 



11 /q-f^Jo) 
M2,S2;M1,S1 V i2 
22 /<-,-(''io) 



.(77,0 ^ 



7 

M2,S2;A»1,S1 



(Ti:-\x2;T::-\x,) 
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The coset decomposition of (j4.1cp needs further, block diagonal self-energy matrices of densities or 'hinge' 
fields (^S^.^^ ^^.^^ ^^(T*2^^\x2;'J'ji^^',a;i) with equivalent structure (j4.2a|4.2bp as the matrices ()4.1al4.1bp in 



order to perform the factorization (j4.3p with the coset matrices Tjjj" ,,2;/ii,si('^i2^^')^2;'J'j/^^ ;3?i) (!4.4ap . The 
generators Y|^2,s2;mi.si('^j2^^''^2;'J'ji^^\^i) (|4.4bp of latter coset matrices consist of two anti-symmetric sub- 



JVi-, ) 



generators Xfj^^^s2;iMi,si{'^ll^'^\x2;7]l^^' ,xi), Xjt2,s2;/ii,si ('Jj2^^\ ^2; 'J'ji^^', ^i) (!4.4cP which describe the anomalous 
paired fermions of a SSB under inclusion of excitonqfj 
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(^-\f2;T;:-\xo)'; 

X2;l' 



(4.2a) 
(4.2b) 



^^/S,s2;mi,^i('^J2'' '^2;'J'ji'' ,^i) = (un-scaled, discrete sums X3,4, jsA = 1, . . . ,2iV, ^3,4 = e, /i, 53,4 =t,i) 



(4.3) 



T* 



i/,s';ii,s\-' j2 iX2i->i^ TXlj 



^ -'At2,S2;Ati,si I ''i4 '-^4, ,X3 



-| ba 
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•(''J2 ) ^ . T-(''Ji ' 
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X;/,s';^,s('J'j2^^ ; ^2; "Jji^^ ; 2?i) ) ; (with complex- valued entries). (4.4c) 



.(1,0 ^ 



(4.4b) 



•^We apply the standard summation convention over double occurrence of indices in (|4.3p with un-scaled sums of discrete, spatial 
vectors £3, X4, and with summations of time indices ^3,^4 = 1, . . . , 2N including summations of electron-, hole-bands /is, /i4 = 'e', 'h' 
and of spin states S3, S4 =t, This kind of summation is also taken within the coset factorization in later steps of transformations. 



24 



4 COSET DECOMPOSITION WITH DENSITY- AND EXCITON-RELATED PARTS 



In analogy to [20] we accomplish the HST (j4.5p by one half with the densities •5^2,s2;mi,si ('^J2^^ '^2;'^'^/^ ,2^1) 

(|4.1ap and by one half with the anomalous parts of dJ^XX-^,, {7^1'^ \x2\ T^'^^ \ f 1 ) (I4.1cl4.1f|4.1gl ) . The 'hinge' 

fields 5S-,,,.^^,,^(T;^'^^',f2;^'^\fi) (!4.1d|4.1eD or <5S2)%„,,;^„,, (T^,'^^', ^2; fi) (|4.2a|4.2bD only have the 
final effect of normalized Gaussian factors (cf. [20]). The true HST for the SSB to coset degrees of freedom is 

therefore given by the integrations of the densities •sJi2,s2;Ati,si ('^i2^^ ' ^2; 'J^,/^ , xi) (|4.1ap and by the integrations 
of the anti-hermitian related 'Nambu' parts of the total self-energy matrix, each contributing in half for the 
transformation of the quartic interaction of Fermi fields (S^*^ = S^"^ q^, cf. ()3.18p ) 
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As we substitute (j4.5p into the original path integral (j2.15p and integrate according to relation p.4p with the 
doubled anti-commuting fields (|2.11al?.llgP , we attain the square root of the determinant and the Gaussian fac- 
tors of the self-energy matrices which are defined in (j4.1all7Tg| ). A shift of the off-diagonal, anomalous-doubled 
self-energy components (j4.6b|4.6cp by the off-diagonal, driving dipole terms 't^^^/.^ ^ (2^2 > 'Jjj^^ ' ; ^i; '^ii''^') (|2.14cp . 

'^^!s';M,s(^2,'J'j2^^^^i,'J'ji''^') (!2.14dp only leaves the block diagonal Hamiltonian parts ()2.14a|2.14bp within the 
determinant and results into a modified Gaussian with a generating dipole term for the quadratic self-energy. 
In order to distinguish between block diagonal 'b = a' and off-diagonal 'b 7^ a' parts, we define in brevity the 
two one-particle operators ^^^?;^.s(^2, 'Jja''^^ ^i, 'J'^i^^') (!4.6dp and ^!!'t>.s(^2, 'J'j2^^'; ^i, 'J'^i"^') (14.6eP for sepa- 
rate density-related and anomalous-related terms of the total one-particle operator H^"s'.^ 3(2^25 'J^j^^'; 2^1, t]^^^') 
(I2.14all2.14dp 
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4.1 SSB with 'hinge' fields and anomalous-doubled parts 
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Ml, 2, Si, 2 2N 



^1.2 il,2 = l 



(''11 ) 



) X 



X TV 

a, 6 



s 

X DET 



M2 ''i2 "'^fl2,S2;fJ.l,Sl\-'j2 '-^2, Jji , J-i; iJ/^i '/jl "^Ml,Si;At2,S2>-''Jl 



V -J ii ; -f- 1 ) -J .; 



1 
A?7 



(z: n-(''j2'. - rr("3l)\ fa 



+ Vj2 'Jv,s'\ii,s\" i2 



02 



X2) 

n '^1) %-i + 



At 



ib=a 



-16a 



T 

J2 



9m ^: 



ba 



1/2 



12 



u,s';iJ.,s\'' j2 ''^2)''ji )-^i;~r 



Qu 1^2 [^AtJ '^iy,s';fi,s{x2,'^j2 ^ j^lj'^ji ^ ; 



ry21 /rr(''i2) ^^.1*"^) ^ \ 



ry21 |'rr(''j2' - 

''^i/,s';/i,sl-Jj2 



^2;T'^\fi)+ 



9m %1 



9m % i 



(4.6b) 
(4.6c) 



(4.6d) 
(4.6e) 



|_|0=a /-f rr'-''J2' . z; rr^ 'Jl'N . 

The coset decomposition follows from the factorization (j4. 3114. 4c]) where one has to consider the integration 
measure from the coset metric tensor and the Vandermonde determinants of eigenvalues [25] in order to replace 



the Euclidean integration variables ^(t]^^^', X2; 'J'ji'^^', a^i)] [El [T71 [20] . A particular simple change of 



integration measure arises from (14. 7j) as one diagonalizes the sub-generators Xf^2,s2;tJ.i,si{'^]2'''^' '^'^''■^h'^ '^iJ 

.„ ^y,^ (T^^- \ X2 ;Tf:«\ XI ) (1133) of coset 



X|(2,S2;M1iS1 V-' j2 



■ M2,S2;M1i'51 ^ J2 



(7-1'^' ,X2;7]l'^' ,xi) (lOcD within the total generator Y'"' 
matrices f |^2,s2;Mi,si('^i2^^'' ^2; 'J'^i^^', ^i) (14.4al) with eigenvalue matrix Y^.^^.^^^(T^ 



■(Vi) 



vectors P"°<,/.^ ..('J'j'a^^', 2^2; xi) ()4.9ati4.9t1) . We regard the anti-symmetry of X^2,S2; 



and of its hermitian conjugate XJi2,s2;Mi>.5i ('^jV^^' ^2; 'J'ji^^') a^i) (I4.4cp by quaternion eigenvalues fc,{7j'^',x) 



n 

x) (j4.8all4.8cl) and eigen- 

Jl 



Mii'^i 

{7^^'^',X2;7]'^'^',xi) 



(S = {+1,-1}) with Pauh matrix {f2)u^ and block diagonal 'eigenvectors' P""'-/i s('^j2^^^ ^2; 'J'^i'^^', a^i) (|4.9al - 



lupt ana diock aiagonai eigenvectors ^Z"s'■^l s^-' 12' 1^2; -^ji"' 
I4.9fp . Furthermore, we use quaternion matrix elements f^^"^^, ^{7^^^^'^\ X2\7\'^_^^'^\ xi) (|4.8cp for the sub-generators 



X 



fT^'^^^ xq-t'"'!' XI ) X'f 

_l-^i2 '-^2, Jj, )XiJ, 



(t'^^-^',x2;t;;^-^\xi) (naEi). 



M2,S2;mi,si >.-'i2 '■^^'•'ii ' -^1^' ^^M2,S2;mi.si ('^i2 ^ iX2]7j^^ ,Xi) ( | 4.4cP of Y|^2,S2;mi,si '^''i2 '"-^'"n 

The correct number of independent field degrees of freedom is maintained by setting the quaternion-diagonal 
elements ^D;u,s;iJ.,s{'^'^i''\^'^ '^j"''' 

{77 } ■ ) 

(T;''2 ,X2;T;/i ,xi) (a = 1,2) (|4.9cl4.9fP identical to zero; these 2x2 quaternion degrees of 



, x) (|4.9ep as part of the total generator ()4.9dp for the 'eigenvector matri- 

ces' P-,. 

freedom alon^ the main diagonal are already contained within the anti-symmetric, quaternion eigenvalues 
'^D;u,s;fi,s{'^j^^ ix) (|4.8bp with Paufi matrix {f2)v 



Y 



ba 



,(7^^\x2;7;;^'\xi) 



y-l;bb (rr■^'^j2^ -TJ^-t'"^' 7\ V^a 
u,s';ti2,si^ 32 T^2i-Jj T^) ^ D;/i2,si;Mi,si 



-paa 

-'"^miiSi;m 



^s{7j , x; 7j_^ jXi) 



(4.7) 



Y* 



ba 

-D;M2,s;Mi.s 



(S^'\x) 



^\);ii2,s;tJ.i,si'^3 ^ ^X) 



^D■,fl2,s■,^il,si7j ^ , x) 



ba 



(4.8a) 
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XD;M2,.;m,^(^^''''.^) = (^2)m2Mi /s(^^\^) ; (S = 2s = {+1, -1} for S =t,i); 



(4.8b) 



fe=0 



(4.8c) 



paa (rr-^^j2 ' ^ . n''''Jl ' z; \ 

p22 /rr("j2) ^^.t'"^!^ ^ 

pll (-rrC^Ja' ^^.T-'^Jl' 

SD;.,s';M,s(^^'^\i'2;T]^\xi) 

S D\v,s;fjL,si^ j TX,7j , x) 
p22 /rrf^w' - .rrt^il) - N 



pll /'n-(''i2) - .rpt^il) -? N 



pll;-l .rj-{''.2) f^.^'"^!' 



p22 |'<y(''j2) - .rrSl) -? X / 

u,s';ij.,s\'' J2 i-^2-i'>ji 1-^1) / ; 



; {11, V = (e)lectron, (h)ole) ; 

exp|-z S5;/.4,«4;/.3,«3 (^J'*' ^4; T^g^'', ^3) } 



i/,s' 



(4.9a) 

(4.9b) 
(4.9c) 

(4.9d) 
(4.9e) 
(4.9f) 



The chosen kind of eigenvalues (l4.7M.8b]l and parametrization with quaternion matrices {fk)v^ (14.8cP takes 
into account the strong Coulomb attraction between electrons and holes resulting into exciton quasi-particles 
of neutral charge. According from a rigorous, mathematical point of view, one could also prefer eigenvalues 
'^D■,^J.,s2■,^l,sl{'^j \ x) = (r2)s2Si //^(•^j ^ j^) with separate eigenvalue parameters /^(T, ^ ,x) for the respective 
electron- '/i = e' and hole-band = h' which are composed of opposite spin pairs due to {t2)s2si for separate 
degrees of freedom of " electront-electronj," and " holet-holej," pairs. This latter kind of parametrization of 
eigenvalues is more appropriate for superconductors with Cooper-pairs of equal kind of fermions caused by 
an attractive interaction, e.g. by phonons. However, corresponding to the strong Coulomb interaction in 
the presented case defined within section [21 the excitons are certainly the prevailing, appropriate kind of 
quasi-particle formation for the eigenvalues XD;^2,s;Ati,s('^j''^\ ^) = {T'2)^i2fj.i /s ('^j''^ ' > ^) with parallel spins of 
electron-hole pairs leading to the two possible, total spin values S = 2s = {+1, —1} for s = {t= +^ , i= — 
One might be astonished about the missing, total spin value S = of electron-hole pairs, having opposite 
spin directions, within our chosen parametrization of eigenvalues; however, the direct generation of electron- 
hole pairs with opposite spins is negligible in comparison to the case with parallel spin directions according 
to the massless, spin-one property of the electromagnetic field with the only two possible spin components 
S = {+1, —1}, but with lacking spin zero component (cf. ()2.3cP ). 

Further factorization of the block diagonal self-energy densities (5S^.^ ^/.^ ^(Tj^^^', X2; t]^^^', xi) (I4.2al4.2bp 
(as 'hinge' fields for a SSB) is specified in relations (j4.10ati4.10fp with the anomalous-doubled eigenvalues 
6A^^A7f'\x) = {+5Xf,47f'\x) ; -5A^,,(t]''^\ x)} (lilObl) and 'rotation' matrices Qif,/;^,,(T,^'', ^2; tJ^^^', xi) 
(liTTna . Qll,.^^^{7^^^''\x2;7^^^''\xi) KWdh so that we regard the hermiticity of St^j^.^^^,.^^^{7^,l''\ X2-,7^^^''\ xi) 
(j4.2ap and their relation ()4.2bp of transposition. The generator 3'd;i^,s';a',s('^j2^^ ;^2;'J'j/^ ,^i) (|4.10cl4.10dp of 
QJJJ"/.^ ^(Tj2^^\ X2; 'J^li''^'; 2?i) consists of a hermitian matrix (I4.10ep . but with vanishing diagonal elements (|4.10fp 
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(77 ■ ) {v ) 

in order to achieve the proper number of independent fields as in <55]^^^, (Tj^^^ ,X2;7^^^^ ,xi) 



(4.10a) 



Aaa;-1 (rr''^]2^ ^ .n-'-^i-i^ ^ \ ( ^"^^3 ,S3 ('^Js ^ '^s) 



a=l 



a=2 



expjz \x4;T]3'3\x3)| 



(4.10b) 



(T]^^^-^\x2;T]:^'^\fi); (4.10c) 



exp 



/^4,S4;M3,S3 V J4 



•J .)4 ' ''^4 ; J J3 ) 2:3 



.('7j2 ) 
i/,s';fi,s 



3"fl;M;.;/.,.(^^\^;^^\^') = ; (;U = c,h; s =t,|) . 



(4.10d) 

(4.10e) 
(4.10f) 



nant for the change of integration measure from Euchdean integration variables d[6T,''J^^, . ^ ^{7]^^^' ,X2;7'j^^'^' ,xi)] 



In terms of the above parameters ()4.7ll4.9f]) and (j4. lOaM. lOf p . it is straightforward to derive the Jacobi determi- 

to the coset part (i[(T~^(iT)^"^,.^ s('^j2^^ > 2:2; T^^^^ ' , xi)] under inclusion of the block diagonal self-energy densities 

(i[(5E^j^ ^/.^ ^(tI^^^', X2; 2^1)] and Vandermonde eigenvalue polynomial Poly(5A^^s(T]''^\ x) ). According to 
appendix A of Ref. [16], we can state the transformation of integration variables 



fj,,u=e 



h {^1,2} 



(4.11a) 



d[<5St».^,,(T]^^«\x2;T;:-',xi)] ^ J] 11 4'JSt>;M,«(^" '^2;^^',^!: 

S,s'=t,iji,2=l,.--,2Af 

d[5±%u,s'-,,A'^^''\x2-,7^''\xi)^ X d[(T-idT)t;^.^,,(T]:^-^\x2;T;:^-^\xi)] ; 

{^} 2 
Poly(5V(T]''^\x)) = n n (5Ae,.(T]''^\x) + (5Ah,.(T]''^\x)) X (4.11b) 

s=t,4.i=l,-,2Ar 

/ {^1,2} 

X I (except : (s' = s) & (ji = J2) & (^i = ^2)) x 

\s,s'=t,iji,2=l,-.,2Af 

X (5Ae,.(T];^i\xi) + 5Ae,..(T;^'^^\x2))' X (5Ah,s(T];^^',xi) + 5Ah,.'(T;^'^2\x2))' X 

\ 1/2 

X (5Ah,.(T];^^\xi) + 5Ae,.'(T;^'^^',X2))' X (ae,.(T];^^',xi) + 5Ah,.'(T;^'^^\x2))' 



4(t-idt)t>;^,,(T,^^\x2;T]:^-^\xi)] 



(4.12) 



g ^d/^!:(T;''^x;7;''^\x)Ad/l;;i,(T;''^x;T;''^\x) ^^sinh(2|7s(T:''^x)|)' 



=Uj=l,-,2Af 



2^ 
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3 /S=2s,S'=2s' {^1,2} 

^ n n n (^p* • =s,s'=s) & 0-1 = & (^i = X2)) 

k=0 \ s,s'=t,l ji,2 = l,-,2Af 
X ■ ■ X 



X 4 



sinh(|/s,a?^\x2)| + |/s(T,^^\fi 



2 I 



I- /^(''j2) - 



sinh(|/s,(<^^',^2)|-|/s(T,^^\xi 



/s.(t;-^',x2)|-|/s(t]^\xi 



1/2 



As we combine the change of integration measure (|4.11all4.12|) with the coset decomposition (j4.2all4.4c|) and its 

eigenvalue parametrization (j4.7ti4.9tp for the total, anomalous-doubled self-energy STP^^,.^ ^(T^^^^', X2; 'J']^^^', xi) 
and also perform the projection onto the coset field degrees of freedom, we finally obtain the path integral 
(|4.13|) replacing (|4.6a|) with Euclidean self-energy variables 



,xi)\ X exp 

O2 '^2,^ii 
-1 ^T-(''J2' ^ .T^^h 'l ^ 



Mi,2,si,2 27V 

if E E v-^(-2,T^:-^x„T^:^^'^ 

^1,2 il,2 = l 



(4.13) 



^j2 9aj2 ^^(2,S2;Mi.siV'^i2 '^2;Jji ' -^iJ ^il ^Ati ^/ii,si;/i2,S2'^ Jl '''^l''^J2 '■^2j 



d[(T-i ^iT);^;,^^_^(T,^^^^^\x2;T,^:^^\xi)] 3 T-i;T;?^ 



X 



det( 



1 

a7^ 



ba 

(• 

v,s'\ii,s 



1/2 



Apart from the Gaussian factor with the density-related self-energy Si]}^,.^ ^{7^!^^'^ ^X2]7^_^^^ the path inte- 



rs "n^ 



gral (|4.13p comprises the functional SIT"-*-; T; %] (|4.15|) with the driving interaction 'Hj'^.^ s(^2, 'J^2^^ ' ^i' '^n 
()2.14cl2.14d|4.6e]) . creating exciton parts, and the square root of the determinant; this determinant has a 
density-related part (I4.14ap . which is completely block diagonal in 'Nambu' indices '6 = a', and a gradient term 
()4.14bp with the coset matrices, which is weighted relatively to the density part (|4.14ap . Furthermore, the source 
term (I4.14cp allows to track the original observables, given as combinations of even-numbered anti-commuting 
fields, to the presented case (j4.13p with even- valued self-energy matrices 

^[^A ^2; fi) = ^^=?;^,,(X2, t]:-'; XI, V (4.14a) 

+ 8''='^ ry,, {i^) 4=?^^,,(T,J^\ ^2; t];'-', xi) q, Vn i 



1 11 



.(I,, ) 



h,5 



ba 

(7^'\x2;'jf;''\x^= (4.14b) 

u,s';fi,s 



rp-l;65' frr^Vj^) .^{1^4) N Ar 



C^Ji '•^4;'J'j3^ 5X3) T„ „ ^„(Tj^^ ,X3;7j^ ,xi)-|- 



ba 



b=a 



|t;;-',x2;^;^^',xi); 



(Jj ,X2,Jj^ '^^ij-l , (J,2 ,X2,J,4 ,X4j 



(4.14c) 



4.2 Saddle point approximation with various kinds of gradient and l/Mx expansions 
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-l:a'a' 



T-i;t;?^ = /45S^i,,,.^,,(T;^-\x2;T;;^^^ Poly(<^A^,,(T^', ) x expj f ^ ^ x (4.15) 

ii,2=i 



2Ar 



X V-\x2y;,''';x^,7l';^^') X Tr 

a, 6 



■^tt^! Via 



/J2 132 /i2,S2;M4,«4 



X 6t 



Ml '^'1 Ati,si;At6,S6^ Ji 



(T;;^-^\fi;T;;;^«\f6) x 



11 . 2Ar 

(T2^«\f6;^^^\f5)T-^^^^^^^^^^ X exp -i Yl 

jl,2 = l 



a, 6 



-'-M2,S2;M4,S4V-^i2 >-^2,-Jj4 ) -^-4; OZj^.^^^^^.^g^ggi^Jj^ ,X4,Jj3 , X3j X 



4.2 Saddle point approximation with various kinds of gradient and 1/A4 expansions 

The factorization into coset and density-related field degrees of freedom with corresponding factorization of the 
integration measure allows to separate the total path integral (j4.13p : corresponding to this coset decomposition, 
the anomalous-doubled fields of the generator ()4.7ll4.9f]) for the coset matrices are regarded to be 'immersed' in 
density-related self-energies of a background functional ()4.17p . This 'background averaging' (I4.17P of the path 
integral (j4.16p for the coset degrees of freedom is denoted by ((. . .))go,a. 



Z[3] 



X DET^ 



(4.16) 



ba \ 1/2 

(t;!-',x2;t;:-\xo 



12 

u,s';ii,s 



2N+1 s=t,l 
X j=0 fi=e,h 



X Tr 

a,b 



'(^^\x;7^^\x) 



X BETU^h,sry^" xexpi-l^ E E E ^^(^2, T^^-'; x^, <'-') x 



a=l,2 si;2=t,l 2N 



^1,2 Mi,2=e,/iii,2=l 



^32 ^f^2 ^^J,2,s2■,^J.l,sl\^' 32 '"^Si jXi) rjj.^ q^-^ 5 



aa 

mii'5i;a'2,s2 31 



(t;:-\xi;t;:-\x,) 



X ( path integral of coset field degrees of freedom 
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^1,2 Mi,2=e,h ji,2=l 



J2 



path integral of coset field degrees of freedomj . 

Moreover, it is possible to derive from the background functional (j4.17p of the self-energy densities a saddle point 
equation (j4.18ap ; in consequence one can approximate the path integral ()4.16p of the anomalous-doubled field 
variables by replacing the density-related field variables sJJJ"/.^ ^('j]^^^ , 2^2; 'J']^^^', xi) with definite, fixed functions 

(s)""^/.^ ^(Tj^^^ , 5^2; 'J'i/^ J xi). The latter saddle point solution results from a similar eigenvalue problem ()4.18bp 
as in section 13.21 but without the off-diagonal term of the driving interaction which has been shifted into the 
functional 3[T~^; T; Ti] (|4.15p . so that the computation of the saddle point equation can be reduced to the '11' 
block of densities 



2N+1 n,s 



11 /T-(''32' ^^.T-f^Jl' 



'j-i ; X3) r]j^ 



-1:11 



3" 



.('7,1 ) 



(x4,T^^';f3,^;?'^V (4.18a) 



(t;J^^',x2;^;^^ ,xi); 

M2,S2;M1i«1 



(4.18b) 



ii=o xi 



X2) 



+ Vj2 («^) (s)i)-,/;^,s(T*2''',X2;T]'''^\xi) r]j^ V'Af;M,s(^]i''\ ^i) = Em i^M-M,s'&jl'^\ 
(eigenstate label M with complex eigenvalue Em) ; 



V'M;M2,^2 (^?'\ ^2) ^M;^i,si (^n''' ' ' ^1) 



V-^{x2,'J^^;x,,'J^-') (s)ii,,,^^^,,^(T;:-',X2;T]:-',fi)/Ar. = 



Mgeigenstates 



Em 



(4.18c) 



The saddle point approximation of ()4.16p with (14.18aM.18cP even leads to a further simplification as one 
considers the relation between the block diagonal propagator (^C[io, (s)] )"■*■'"" and the self-energy density (5)"'^ 
as consequence of (l4.18aM.18cP 



-l:b=a 



(4.19) 

After insertion of the saddle point solution (s) from (j4.18all4.18cp with inclusion of (|4.19p . we can remove 
the block diagonal propagator (3i[.^, (s)])^"'"'"'^ of purely density-related terms and obtain the path integral 
(|4.20a|4.20bp with the gradient term (|4.2Ucl4.20dp which simplifies the gradient expansion of the trace-log 
functional. In this manner one has succeeded into a path integral of exciton related coset matrices 



Z[d] 



2N+1 s=t,i 

4(T-^^T):;.^_^(T;;'-\f2;T^\f,)] 3T-;f;7i xexp^^A^EE E >< 

X j=0 ii=e,h 



(4.20a) 



4.2 Saddle point approximation with various kinds of gradient and l/Mx expansions 
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X Tr 

a,b 



Z[3] ^ / 4(t-^ dT)^;,^_^(T^\f,;T^\f,)] 3 



2N+1 s=t,l 



M^-v^EE E S 

X j=0 ii=e,h 



J 5 -^5 "J 



In 1+ 



(4.20b) 



-1 



exp{[Y, ...]_; [^,(5)]] 



5^^^[^,(5);t-i,t] « (exp{[Y, ...]_}^C[^,(s)]) (v-i (s)- /AA,^^ - 1 



(4.20c) 



(4.20d) 



There are various possibilities to expand the determinants in (j4.20a|4.20bp . One possibihty refers to the param- 
eter Mx = {L/Axy of total space points which allows to consider the self-energy (•s)^'^/.^ ^(^jj^^ ■iX2'i'^ji^ i-^i) 
as a background density provided that its spatial variation is negligible compared to that of the coset matri- 
ces. A further possibility follows from performing the total trace for the determinants in (j4.20a|4.20bp in a 
momentum basis of plane waves with a gradient expansion of spatial derivative or momentum operators [32j ; in 
this manner one can derive effective Lagrangians of coset matrices with a minimum number of derivatives and 
effective coupling coefficients consisting of the saddle point solution. One of the guiding principles for such a 
gradient expansion has to include 'Derrick's theorem' [33] which restricts the expansion to stable configurations 
under a minimum number of derivatives of coset matrices. Moreover, one can introduce a generalized gradient 
operator in the two following kinds 



9^ (matrix)^" : 



^ln(-5f) (matrix 



ba 



^([S), (s)]'' (matrix)''^ (s)]"''""^ - (matrix)''^ 

exp |[ln(-3<[^,(s)]) v]^} - (matrix)^'"' 



(4.21) 
(4.22) 



where the exponentiation — ^[i^, (s)] = exp{ln(— ^[i^, (s)])} allows for a further possible expansion in terms 
of [ln(— CK[io, (s)]) , . . .]_. Since the one-particle operators -^""/.^ ^(2^2) 'J'jj^^ S^ij'^ii''^ ) contain the total unit 
operator (in fact '—1', cf. (j2.14al2.14bp . therefore d^^^_^^ instead of d^^^^^ in (j4.22p ). one also achieves reasonable 



ln(3<) 



approximations and simplifications from the expansion of the logarithm by using ln(l — x) 



1^2 1^3 



ln(-:K[:^,(s)]) =ln( - S)iy^^Jx2,7ll^^'-,x,y;;^'') - q^rj,, {^f) (s)J,^,.^,,(T;;'-\ xs; f i) r?,, 



Jvi-, ) 



In 
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5 TRANSFORMATION TO EUCLIDEAN PATH INTEGRATION VARIABLES 



However, the appropriate expansion or approximation of the determinant in (j4.20a|4.20bp relies on the prevailing 
phenomenon which may refer to universal fluctuations of a 'Random Matrix Theory', an effective expansion 
with second and fourth order gradients for Skyrme-like Lagrangians or extraction of relevant terms for a 
renormalizable theory. Therefore, 'ultimately' valid expressions from expansions of the determinant are hardly 
accessible or universally applicable so that the prevailing approximation of terms has to be adapted to the 
corresponding physical phenomenon which has originally been defined as a coherent state path integral ()2.9al 
I2.10cp [6j, a Dyson equation or further diagrammatic expressions [7]-[l5]. 



5 Transformation to Euclidean path integration variables 
5.1 Transformation with the square root of the coset metric tensor 

In the following we abbreviate the contour time variable T^^''* and the spatial vector Xj by 'J'^/^. and also 
combine the total set of variables with electron-hole and spin indices into respective numbers 'i', '2', . . . for 
brevity 

'^j/' ' = "^ji^i ; (5.1a) 
'1' ; '2' ; . . . ^ '(t];^'^' , XI , , siY ; '(t]^^^'^' , 5^2 , /X2 , ^2)' ; • • • • (5.1b) 

Section 11.21 contains a description of a nontrivial, multi-dimensional integral with the square root of a metric 
tensor as the Jacobi matrix. In order to convey these transformations to coset spaces as S0(9T, 9T) /U(9T), we 
identify the matrix combination T~^(aT) as the analogous square root of the metric tensor as in section [L2l 
The derivative symbol 'a' can denote a partial derivative of space ^d^ = ^' or time 'S/T*''^'', a total derivative 
'd' with respect to time slices for the contour time ordering of path integration variables or a variation symbol 
'5' for deriving classical equations in the exponentials on the condition that the integration measure has been 
transformed to Euclidean variables as in this section. The derivative symbol 'a' can even refer to the matrix 
operations (|4. 21114. 23|) for generalized derivatives (dj^T), (c^in(_j{)T) of the coset matrix (j4.4all4.4cl) . We achieve 

the transformation to Euclidean variables by using the property that 'T~^(aT)' can be expressed as elements 
of the so(9T, 9T) Lie algebra. This can be illustrated by consideration of a Lie group U ((/?*) of matrices with Lie 
algebra generators ti and parameters As we take the variation AU((p^) = U{ip^ + A(^*) — U{ip^) of U{ip^), we 
can apply the closed Lie algebra of generators ti in order to relate U ((/?* + A(/3*) to a rotation of U ((/?') by another 
Lie group element U{Af^{ip^,Aip^)) with function increments a/-' (93*, Ac/?*) ) to be determined. Nevertheless, 
we can specify C/^^((^*) {AlJ^ip^) ) in terms of the generators tj 

U{ip') =exp{iip' ii} ; (5.2a) 

U{ip' + Aip') = exp{« {ip' + Aip') ii} = U{ip') U{Af{ip\Aip')) = exp{i ip' ii} exp{i Af {ip\ Aip') ij} ; 

(5.2b) 

U-\p^) (AUip^) ) = U-\p) {U{^' + A^') - U{p)) = t Afi^\ A^') i, = I a/ ij . (5.2c) 

Application of (|5.2aM5.2"c]) to T~^ (aT) combined with an additional gauge transformation P . . . therefore 
gives rise to the so(9T, 9T) Lie algebra elements aZ''"('5'; '1') (I5.3ap with the block diagonal, density-related 
parts Ay^^('5'; '1'), Ay^^('5'; '1') (|5.3d|5.3ep as dependent variables and off-diagonal, anomalous-doubled parts 
a!X('5'; '1'), aDC^('5'; '1') (I5.3b|5.3"c]l as the independent variables which we also define by quaternion-valued ma- 
trix elements ^f^issi ('J'jg''^g j 'J'j^''^|)- The block diagonal density parts consist of anti-hermitian sub-matrices 
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I Ag('5';'i') which are functions of the coset variables aX('5';'i'), aX^('5'; 'i') or their quaternion-valued corre- 
spondents ^^.(7^^;'!';^]^^, ^^(^^^^^^^^^ 



_ / Ayii('5';'i') aX('5';'i') 



ba 



P('5';'4') T"^('4';'3') (At('3';'2')) P'^{'2';'i') 

3 

Aa('5'; '1') = y (ffc) Aai^) (T^^V ; T 



aX('5';'i') — aX^5_s5;^i_si('J^.^^|_.;'J^.^'^1^ 



fc=0 



Aa('5';'i') = -Aa^('5';'i') ; 
Ay^^('5'; '1') = -aV^^'^('5'; '1') = z A5('5'; 'i' 
A5('5';'i') = A5t('5';'i') . 



(5.3a) 

(5.3b) 

(5.3c) 
(5.3d) 
(5.3e) 



In order to calculate the dependence of Aa('5'; 'i') on the anomalous-doubled fields within the generator Y('5'; 'i') 
(|4.7ti4.9fp of coset matrices T (|4.4ati4.4c]) . we use the property 



exp{i?} (5( exp{ - B] 



dv exp [v B] 5B exp{ -V 13} 



(5.4) 



and perform the derivative symbol 'a' of P T ^ (aT) P ^ with its additional gauge transformation so that the 
eigenvalue decomposition of Y('4'; 'i') into P^^('4'; '3') Y£)('3'; '2') P('2'; 'i') simplifies the expression considerably 



a2,^''('5';'1') = (P('5';'4') t"^('4';'3') (At('3'; '2')) P"^('2';'i' 



ba 



(5.5) 



P('5';'4') exp{Y('4';'3')} (a exp{-Y('3'; '2')}) P ^('2';'i') 



ba 



ba 



P('5';'4') (exp{P"^ Yd P}('4';'3')) (Aexp{-P"^ Yd P}('3';'2')) P"^('2';'i'; 



dv (exp{?j Yd}('5';'4')) (aY'('4'; '3')) [exp{-v Yoji'^';'!')) 



Note, however, that we do not directly transform from (aY('5'; 'i') ) ()5.6ap to a2,^"('5'; 'i') in (j5.5p . but from 
rotated variables (aY'('4'; 'i') ) = P('4';'3') (aY('3'; '2') ) P*^('2';'i') (I5.6bp to a£''"('5'; 'i') which does not alter 
the integration measure due to the invariance of the trace (|5.6cp 



aY('4';'i') 



aX('4';'i') 
AXt('4';'i') 



^ 6a 3 

) ; AX('4';'i') = y(rfc) a/^^^ (t'^^V ; T^.^^V ) ; 

^ k=0 



(5.6a) 



aY'('4';'i') = P('4V3') aY('3';'2') P-\'2';'i') ; A±' {'4' ■'!') = ^^^^ ^f'tS.sS^&'^h^^ ^ (5-6b) 



fe=0 



Tr 



aY'('3';'2')) (aY'('2';'i')) = Tr (aY('3';'2')) (aY('2';'i')) 



(5.6c) 



It remains to compute the exponential (exp{zbt> Y£)}('5'; '4') ) of the quaternion-eigenvalue matrix Y£)('3';'2') 
(14 .7114 -Of]) with anti-symmetric element {t2)u^ and to integrate over v G [0, 1] in (I5.5p in order to relate the 



34 



5 TRANSFORMATION TO EUCLIDEAN PATH INTEGRATION VARIABLES 



'rotated' independent fields (aY'('3'; '2') ) (lOT]) or a/^SJ^^^CT^^'V ; T^^'V ) to aZ'"^('5'; '1'). This has to be ac- 

comphshed for the various parts a, 6 = 1, 2 of aZ^"('5'; '1') ()5.3all5.3^ where the off-diagonal parts comprise the 

independent integration variables Aaf^^^s5;iii,sii'^j^%^','^j^^g^) (|5.3b|5.3cp and where the dependent, block diagonal 
parts are to be determined by corresponding off-diagonal entries. We emphasize that the eigenvalue decom- 
position Y('4';'i') = P~^('4';'3') Yd('3';'2') P('2';'i') (j4.7ll4.9tp diagonalizes the coset metric tensor in analogy 
to section [L2] in a 'transferred sense' and yields straightforward, simplified expressions. Therefore, we simply 



(|5.3bl5.3cp in terms of a/'^''^ 



J (I5.6ap and compute the respec- 
tive change of integration measure so that we can attain Euclidean integration variables Aoi^s^ '^'''^ ' 



^5,3:5' jl,Xi' 



list the various parts, following from (15. 5p . and specify the various variables Aa('5';'i'), Aaig^ 

)n 



^5' jl,xi'' 



under consideration of the integration measure (14.121) for Af^.i „ :7^^'i ) or AfLl, , (T^^ 

diagonal elements for the anomalous-doubled (15.3b|5.3"c]) and block-diagonal terms (I5.3d|5.3ep are given in 
()5.7ati5.7"c|) and ()5.9p with the respective change of integration measure ()5.8p which cancels the corresponding 



term in ()4.12p and results in Euclidean path integration fields dosl'' (Tj V ; t] V ) A dci~ss* ^^'-'l' 



P T-^ (aT] P 



-1 



P T-^ (at) P- 



(2), 

■ss 

21,t 



\ j.S ' j.x J 



Aa 



SS \ 3,3 ^ ^j,x ) 



^2, 



+ 



sinh(2|/s(T];'|')|) 



(5.7a) 



+ 



1 sinh(2|/s(T;;'|')|) 



4i7s(^:i')i 



exp{.2 0s(T;y)}A/2;;(T<^';Tj::£') 



( 



sinh (2|/s(t]J)|) 

4|7s(t!1')I 



V 



sinh (2|/s(Tg')|) 

4i7s(T*;'i')i 



?/(2) /rr(';j).rr(';j)N 

?/(2)*/^(^j). rr(';j)^ 
^JD;ss ' •'j.S J 

/ 



I7s(^)l 



sinh 2 



i7s(^::^')i) 



sinh (2 



=^2<^s(Ti::i') 



sinh (2|/s(t';'|')| 



4 i7s(^;i' 

1 I sinh(2|7s(Tj;'i^|) 

2 4|7s(t];'|')i 



sinh (2 |7s(t':'|')|) 

i7s(t;:^'i 



sinh(2|/s(T];'|')|) 



(5.7b) 



?/(2) 
p'(2)*/ 



^•1 Diss 3,3 ) •'j.S ^ 



(5.7c) 



Aai?(J^y;T;:i') 



Jf(2) /rr(''j).rr(''j)N . J f (2)* /rr(''jV. rr^"j 
"■J D;ss\-' i,3 ' •^j.S ^ "-J D;ss\-' j,S ' -^j.S 



r('7,) rr('^^)^ 



2£s(Tg) 



sinh 2 



l7s(Tlf)l) 



(5.8) 



P T^ 



AT) P- 



V i,2 ' 3,3 ) 
/X5,s;/ii,s 



22 T 

Pt-ifATlp-M ' (Tfj^,ti) = ^^9,,,s,,usi7tS^'^t^)= (5-9) 

/X5,s;/Ji,s 
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:sinh(|7s(T;g'|)) 



-6 i 



tanh(|/s(T3')|) 



Similarly, we calculate the off-diagonal matrix elements for anomalous-doubled (|5.3bl5.3cp and block diagonal 
parts (j5.3d|5.3ej) . Expressions are abbreviated by coefficients A{'5';'i'), B{'5';'i'), C('5';'i'), D{'5';'i') in order 
to attain corresponding results from calculating and analyzing matrix terms in eqs. (|5.3ati5.3e|5.5l5.6all5.6cp . 
Note that we again obtain the result of eqs. (j5.7all5.9"]) with the diagonal elements, as one takes the limit process 

('7 ■ ) ('7 ■ ) 

'5' '1' of the spacetime labels '^'j^^lg ""^jihi '^ith identical spins S5 = si and 2x2 quaternion elements 



P (AT) P-)" (S^i^l^^J -(P t- (At) P-)"'^ C^i^^:^^ = (5.10a) 



^('5';'!') 



|7s.(^S)hosh(i7s,(a;::i;)i).nh(i7s,(T;3^)i)- 




cosh(i7s3(T;3;)i) sinh(i7s,(Tj;fi;)i) 


|7s.(^:i;)r- 

|7s5(^:5;)i cosh(i7s,(^:i;)i) sinh(i7s,(T;::i;)i)- 




2 

cosh (|7si(Tft^l ) sinh (|7s5(^:i^l ) 


l7 fT^'^sM'- 




2 



(5.10b) 

l7. rT<''^5M _v,ri7. rT<''>5\n «inhri7. rT^^^iNiUlx tt'-'^iM .n«h ri7. (T^^n^A ri7. rT'''^5\,^ 

5('5';'i') = 



Ao^^i,(TV'^l';T*.''V^ 

J5,X5' ji,a;i' 

Aai°l:(T*.''^i^T*.''^V^ 

SbSiK j5,a:5' 



(5.10c) 
(5.10d) 



ACsvi') exp{z(^s3(T;::i:) +0s.(t;::i;))} i^csvio \ / ^m.s.&;:^'^;& 



V expf z(0s3(^:i) +<^s.(a^::^:))} vio acsvio ; ^ 

^('5';'!') -exp{^(0s3(^:i;) +0s.(t;::i;))} i?('5ri') \ / Aai^k(^:i:;T;^,,^ 




expf ^(0s.(T;::l;) +^s.(t;::i;))} b['s'-'v) Ai'sri') J V 

(5.10f) 



Aaitk(aj;;i;;T;:^i;) \ .=1^2,3 



AaiS:(T^.''^l';T^.''^l' 

/ ^('sri') -exp{^(0s.(T;::i;) +^s.(^:v))} ^^csri') \ / A/s^uiTt^i^'j^^: 
[ -expf ^(0s.(^::l;) +^s.(t;::i:))} scsvio acsvio ) [ ^fSL^&'^^:^ 



(5.10g) 
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7 (t'^^S 


+ 


7 fT*"^!') 




sinh 




+i7s,(^::i;)i) 



^2e5vi')-52('5';'i') Hsinh(i7s,(<i)i-i7s,(^:^;)i) sinh(i7s,(a^:i)i+i7s,(^:i;)i)y; ^' ^ 



— rlnC^) CT^''J5' . rr'-^'hU . J-(fc)* /rrt^JS^ . rr'^Jl ' \ l7s5 CJ^jg^ ) l + l7si C^'ji^l^ ) I l7s5 (^js'^lg ) I" l7si ("^0^% ) I . 



(P (At) P-)" iT^i^-y^i) (P (At) p-f (7^i;7^i) A^Csri' 



k=l,2,3 

C('5';'i') = 



(v ■ ) 

^^5^ J5,a:5^ D('5';'l') A f n ^ (T- - iT- - 



(5.12a) 

+ 



V ' / •'D;S5SiV J5,a;5' Jlja;!' V ' ^ ■) D;S5Si^ J5,xs^ Jl,xi' 



(777,2) = (to) ; (ihs) = —I (fi) 



(5.12b) 
(5.12c) 









) cosh(|/sg(T]3;)|)- 


/S5(^S.!)| (i/si(T];:i;)i) sinh(i/s,(Tji;)i) 









i:)('5';'i') = C('i';'5' 



(5.12d) 





7s5(^i)k 


7s5(^S)i-h(i7s,(a^:i)i) cosh 




- 


7s.(^::i;)i sinh(i7s,(c/:^^)i) sinh(i7s,{cr;:«;)i) 







11 



pt-i(at)p-M (t; 



P T"^ (aT) P 



-1 



M5,S5;/*1,S1 



(5.12e) 



-«h(i7s,(^:S)i)+cosh(i7s,(^::i;)i) 



+ 



cosh (l/3^(T]3;)|)+cosh (|/s,(T;;^V|) 



(I/s^t'^^^I) AaW•(T"'^5^T"'^l') 

\'''^5^ J5>a:5'v 'S'l^ ^5,2:5' 31,^1' 



5.2 Eigenvalues of coset generators and their transformed, Euclidean correspondents 
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fc=l,2,3 



-J (/>S-, (T, ) ^CO) /o-(''j'-,' rr(''jl K , 

sssiv ^5,2^5' s^siK ^5,^5' ji.xi 



X e 



fc=l,2,3 
i </>s 



— tanh 



+ 



fc=0 

_ Aai'^); (jC.s.) ^(".ij )1 _ tanh ( ^^^^^^^f^^^'^^^^ ) 



X 



e- ''is;' A«iS, ; 3^:5: ) + (-1)' e' "IS' a;.*); ^.^^^ 



} 



+ 



} 



■'''Jr\ , . , ) (J7 ■ ) 1 ~1 -I- 

\ (matrix)' = (matrix) 



5.2 Eigenvalues of coset generators and their transformed, Euclidean correspondents 

In the previous section [5?T] we have accomphshed various relations (j5.7all5.8|5.10all5.1ip between the rotated vari- 
ables (aY'('4'; '1') ) = P('4';'3') (aY('3';'2') ) P~^('2';'i') or their quaternion representation with aX'('4';'i') = 

ELo('^fc)M4Mi^/Dfi4Si('^jI^j'^ji'5;!) and the Euclidean path integration variables l^a^s^sx C^^ ^, ' '^hA [ ^^^'^^ 
from the anomalous-doubled fields, we have also transformed the block diagonal density parts (|5.9|5.12all5.12ep 
as dependent variables on the off-diagonal terms. The eigenvalue terms Y£)('3';'2') or fsC-^j^s ) of Y('4';'i') 
(|4.7ti4.9f]) have diagonalized the coset metric tensor and have simplified the resulting transformations between 
AaissiCJ^]!''!^; T^.^^l^) and A/^'^g^gg^('J^.^^||^;T^.^^l|); therefore, it remains to determine the exciton-related eigen- 
values /sC^j^l') (I4.8b|4.8cp with absolute value \fs{'-^'j\s)\ and phase (/>s(TjV) in terms of the Euclidean fields 
(iai5Si(T^-^^|^;'T^-^^l|). Since the coherent state path integral has been derived from a contour time ordering of 
the time development operator, one has to apply the infinitesimal integration increment 'd' of the contour time 
in order to relate the Euclidean fields da('5';'i') to the eigenvalues df^{7j j; ) 

ifs(m - i7s(7;3)i = xS:'"' ■• *=<^'"'') - *s(TS) = xS;>" ; (5.13) 

a fT'"^'- t'"^^) - a fT(+^- T^+^l -5 ' dT*"^'' ^^'^'"''^ ^'-^ ' '''' ^ (5 14) 

As we reconsider the relation ()5.5p for 'A'='d' of the contour time and expand dY = {dP~^ Yjj P), we note 
the additional commutator [Y^i , ((iP)P^^]_ whose resulting diagonal quaternion elements in the off-diagonal 
parts can be chosen to vanish due to the gauge invariance of the integration measure with respect to P . . . P^^ 
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(cf. (j5.6cp . compare appendix [A|) 



(d£''*('5';'i')) = (P {(if) p-^y^B';'!') = - f^dv (e^^« P {dP-^ Yd p) P~^ e"'' (V; 'i') 



dv 



,v Yd 



{dYo) + Yd , (dP) P 



-1 



:,-v Yd 



ab 



('5';'!') ; 







dP)p-')(7^l^;7^i'); ( Yz, , (dP) P^^ )CJ'^i' ^'^T^^ 







(5.15a) 
(5.15b) 

The chosen gauge (I5.15bj) simplifies (I5.15ap and leads to the direct transformation (I5.16ap between the diagonal 
elements of the Euclidean path integration fields dZ°'^{'5'; 'i') and the eigenvalues dYD where the integration of 
V G [0, 1] and the calculation of exp{±?; Y^)} becomes straightforward due to the restriction to diagonal matrix 
elements 



dZ 



ab 



\''i,s 1 •'j.s J 



ab 



dv (exp{vYD} (dYD) exp{-vYD}y(7^i\7^i'' 

\ / ii,s;u,s 



dZ 
dZ 



12 (cj(^iK <j(^j^ -^ 



21 



A'' 3,3 1 •'3,3 



(5.16a) 

(5.16b) 
(5.16c) 
(5.16d) 
(5.16e) 

After computation of (I5.16ap . we achieve for the various diagonal elements of (I5.16bll5.16ep following relations 
= d{\as(7^i')\ exp{z «s(Tj;^')}) = (5.17a) 



"'^ti,s;i/,s\-' 3,3 ) •'3,3 ) 
JO- 22 



K'' 3,3 ' 3,3 I 



^ dg^^s;v,si!^ '^'^3,3 ) > 

~1' ^fiu dg^^s;v,s^ j,x 'i'^j,3 ) 



dfl(T^^) eM^'i(t>s{t,3)} 



1 sinh(2|7s(T];'|')|) 

^ 4|7s(t];'|')I 



+dfsCJ^i)[h + 



Ida f'T'''^'.'T'"^h 



ytM,s;ii,s\-> 3,3 1 •'3,3 J t 



{dfs(7^^)) exp{-z 0s(T:::£')} - {dfl(7^l')) 4>s(7^i')} 



1 , sinh (2 |/s(T;"|')|) • 
4|7s(^S')l . 

(sinh(|73(Tj-^^|))^ 



2 I /S (^''1^1 



(5.17b) 



and da<^(7^i>;7'^i>) 



da (j''^j'>.'j^^j'>] 
'^yfi,s\^,s\-' 3,3 ) ^3,3 ) 



d{\asm')\ expO as(7^J)}) = 



sinh(|/s(T;y) 



d<^s(7iy) 



(5.18a) 



(5.18b) 



where we have also introduced the separation of the exciton-r elated, diagonal elements dai^^ (t] V ; V ) = 
d(|os('Tj^j^)| exp{? asC^j^j')}) into absolute value |as('J']^i')| and phase as('J'j|2')- However, in order to compare 
real and imaginary parts of the transformations (I5.17ati5.18bp , we perform a phase rotation from da^^J ('J'^'s ^ ; '^'j^s ' ) 



to (iais''('jj|'|'; 'jji'j') with absolute value |as('T]'2^)| and corresponding phase as('J']^'| 



~(2)/rp(')i).rr('7j)^ 



i5s(^::^^) 



exp{z as(x5^)} 



(5.19) 



5.3 Nonlinear sigma model with Euclidean coset fields S0{% 01) / U{m) (g) [/(OT) 
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^ da^g^ > ) ''^ d^ss i'^j.x ^'^j.x ) dcL^J ('J'j £ ) '^j.s ) ^ dalf.^{7j ^ s ) • 



Accordingly, we can replace the diagonal, Euclidean, pair condensate integration variables of eqs. (jS.lTal - 
[5l8b]l by (£i^.^(Tji^Tji'), (+C.C.) and obt ain new relations between the coset eigenvalues fsi'-^j,s ) = 

I7s(^jy)l exp{^0s(^3')} and the new Euclidean elements T^^j') = |Ss(T];i')| exp{z5s(Tjy )} ([EIS]). 

Note, that the absolute values transform as state variables whereas the phase transformations are path depen- 
dent and , therefore, can only result from contour time integrals 



i/s(^si) 



|5s(T}:i')|-|Ss(7J5)|; 



(5.20a) 
(5.20b) 



=0 



0s(t;:i')-<as(^5) 



( + ) / Q^p') 



3 ,x 



3 



(5.20c) 



3 



(5.20d) 



Consequently, we have the final relations between absolute values (|5.21ap and the phases ()5.21bp where the 
Euclidean field variables ai? (T^V; T^V) = |as(T]V)| expjzasCjjV)} have to replace the original diagonal 
elements t]V) by subsequent relation ()5.21cp . due to the inclusion of the phase rotation (|5.19p 



i/s(^:i')i 



^s(t}:i') 



|Ss(^y)l ; 

exp{<</.s(^?) -5s(t5:^'))} d^^is^^'^'Jt:^) 



dT/ 

■{+) J 



i 



sinh(2|S:s(T7 )l) 97,' 

^ 3 ,S ' 3 



exp < I 



dT/ 

•{+) J 





(•I ■/) 

2 |5s(n^ )l 
J ,x 



3 ,x 

3 



(5.21a) 
(5.21b) 
(5.21c) 



5.3 Nonlinear sigma model with Euclidean coset fields 80(01, DT) ® \J{% 

It remains to collect and summarize the transformation to Euclidean field degrees of freedom and their relation 
to the eigenvalues of the coset generators following from sections 15. H 15. 2[ According to these transformations, 
the final nonlinear sigma model, following from (|4.2()all4.2()d]l is determined to be 



X exp 



2N+1 s=t,4, 
X j=0 fj,=e,h 



(5.22) 



X Tr 

a,b 



In ( i - (5(1^ -JT)^)'''"' + P a[t-i(5), t(a)] (s 



b=a 



it 



3 1 -^I 3 1 
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with the sub-functional 3[T ^(a);T ^{a);'H] (|4.15p containing the driving laser field for creation of exciton 
quasi-particles in the off-diagonal, anomalous parts 



3^ 

X V"^(x2,1j2;xi,Tji) X Tr 



t-i;T;7i = I d[6t]^.^,^,.^j7';^^^\x2;7^^^^ Poly(<^A^,,(Tf^\ x) ) x exp^ l ^ J] x 

ii,2=i 



2N 



(5.23) 



a,b 



JVj„ ) ^ 



X T 



-l:a'a 



2Af 



X exp<;-i ^ 

jl,2=l 



X V-\x2,7,,;xi,7n) x Tr 

a, 6 
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M2,S2;M4,S4V" J2 



{Jj ,X2,Jj. 1X4) Ol^D g g {J j ,X4,Jj ,X3) X 



The increments of Euclidean path integration fields Aa('5';'i') (|5.3all5.3e|) are specified by quaternion matrix 
elements Aai^L ('J'!''^^' ; 'J'!''''^' ) with the exception of the diagonal elements Aais^T-l'; T-''^^) which incorporate 
a contour time dependent phase correction (j5.19p . denoted by the tilde ' ' in Aa('5'; '1') 



a5('5';'i') = -A^ ('5';'!') = e(^)msmi 



fc=0 



2|Ss(t;:^') 



\smh{2\as{1pj)\) 
(last term (tJ;:^' := 7^^^ = 7^;^%)^^ := s, = s, , S = 2s)) ; 



3 



4?ro';7^J) = |Ss(^T)l exp{.5s(T%0}; 



4sS(t::^;t^-')] = n n 

s=Uj=l,...,2Af 



S=2s {x} 



2i 



(5.24a) 

(5.24b) 
(5.24c) 



3 /S=2s,S'=2s' {^1,2} 

^ n n n (^^^^p* • =s,s'=s) & (ji = ^2) & (xi = X2)) 

k=0 \ s,s'=t,i ii,2=l,...,2Ar 



,.{fc)*,^{lj2'.rr(''jl)N . J"(fc)/'<y(''j2).rr(''jl)N,\ 1/2 

2 1 I 



Aside from the source field 3[T ^(a),T ^(a)] for generating the observables, 



6a 



i ,i (J ,X2,J,^ ^^V - ^,y,s';f^i,sMj2 '^2,J,4 , ^4) X 



(5.25) 



1:6'6' <i6'a' 



5.3 Nonlinear sigma model with Euclidean coset fields S0{% % / U{% (g) U{% 
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we have introduced generalized gradient operators (j4.21|4.22p . ()5. 26115 ■2"8]) with the adjoint action of the loga- 
rithm of the density-related one-particle operator (|4.23|) 



drfi (matrix) 



ba 



(s)]''^ (matrix)''" (s)] "'"'"'^ — (matrix)*"^ 



^ln{-3f) 



^in(-5£)^('3';'2'))('5';'i'): 



ba 



exp|[ln(-tf{[^,(s)]) ,...]_} -ij (matrix)*'"' 

n^,m ,■■■][} - i)z''''^'('3';'2' 



exp 



ba 



ba 

) ('5';'!') 



(5.26) 
(5.27) 
(5.28) 



The logarithm of the density-related, block diagonal operator !K[S), (s)] allows to simplify the action of the 
adjoint operator adj^^j- in ()5.27l5.28p by expanding the logarithm around the unit value with the saddle 

point solution (s)""/.^ ^('J'jj^l^' '^j/^i) ^^"^ kinetic energy operator of the one-particle part 



In 



J»7,„ ) r^iVi^ ) 



An /i\ll /^(^ul ^ ^ 



11 

/i4,S4;/i3,S3 



-1:11 



2N+1 ^J.,s 



M2,S2;Ati,si 



ii=o XI 



+ Vj2 Qu {ijf) (s)"^';^,^(T> ,X2;7^l'' ,xi) r/j, g^J'0M;M,s(T,i'' ,xi) = Em iJM;u,s'(7^l'^ ,X2) ; 
Mj^h,m) (Tl^^"\x2;T,*:-\fi) = V-i(f2,T,,;fi,T,J (s)^=?^^,,(t;,^-\ x^; T.^', xi) 



(5.29a) 
(5.29b) 

(5.29c) 
(5.29d) 



The coset matrix 2,^'^('5'; '1') (|5.3all5.3e|) . (|5.3Uall5.3Uc|) consists of the anomalous-doubled, exciton-related, anti- 
symmetric matrices Aa('5';'i'), Ao ('5';'i') as the independent field degrees of freedom which prescribe the 
dependent, block diagonal, density-related parts (j5.30d|5.30ep 



a2;'"'('5';'i') = -('P('5';'4') t"^('4';'3') (At('3'; '2')) P~^('2';'i' 



ba 



Ayii('5';'i') Ai('5';'i') 
AXt('5';'i') a922('5';'i' 



ba 



(5.30a) 



AX('5ri') = ax^,,,,;^,,,,(t;.;;i;;t;.;^i;) = M'^rr) = Y.i^k) 



k=0 



Aa('5';'i') = -Aa ('5';'i') ; 
Ay^^('5';'i') = -Ay^^''^('5';'i') = I A^('5';'i') ; Ag('5';'i') = A5^('5';'i' 



(5.30b) 

(5.30c) 
(5.30d) 



42 



6 SUMMARY AND CONCLUSION 



tanh 



k=0 

2 



tanh 



|as.(T*'''5^|-|Ss-,{^*'''l^l 



X exp 



r(+) ^ 





-1/ 



+ 



X exp O / aT / 



(5.30e) 



6s,s, tanh(|5sTO0i) |Ss(t;;:2')I a5s(t;;:£\ , 

(last term (T^^ := tJ.^I^ = := 55 = si , S = 2s) ) ; (matrix)t = (matrix)*'^ ; 

(cf. gauge fixing in eqs. (15.15a|5.15bp and appendix [A]) . 



The transformations to EucUdean fields yields corrections for the derivation of classical field equations by first 
order variation of the anomalous field parts Aa('5';'i') and becomes even more important for higher order 
variation of fiuctuation terms (cf. section ll.2p . The derivations to Euclidean path integration fields within 
sections lS-H 15.21 therefore hint at the importance of the correct coset integration measure, orginally included in 
the factorization of the total self-energy ^S^^*^^,.^ s ('^jl^l2 ' '^i?^! block diagonal, density-related self-energies 
as hinge fields of a SSB and into anomalous-related, off-diagonal parts by using a coset decomposition. 



6 Summary and conclusion 

6.1 Observable quantities by differentiating with the generating source term 

In sections 13. II to [3^ we have performed the three possible Gaussian transformations of increasing complexity 
in order to convert the quartic interaction of Fermi fields to even-valued self-energy variables and matrices. 
These path integration fields, following from the prevailing HST (|3.1ati3.3|) . (|3.9ati3.12c|) . (|3. 16113. 2T]) of sections 
l3.1l to [3r3l only take values in completely Euclidean spaces and flat integration measures so that the calculation 
of observables results from differentiation of the corresponding HST-transformed generating function (13. Sal - 

[33c]l . (I3.13all3.14p . (|3.22ati3.22bl) with respect to the source field Sty.^^si'^py,^^^) dHOaD. The latter source 
field allows to track the original observables, as combinations of anti-commuting fields in even number, to 
corresponding Green functions with the prevailing self-energy variable or matrix which can be restricted by the 
computation of a saddle point equation ()3.6ll3.8p . (|3.15all3.15cp . (|3.23cll3.23fp . 

We outline the calculation for the more profound SSB case with the factorization into density- and anomalous- 
related parts by a coset decomposition according to sections [H and [3 The density- and Nambu-related parts 
follow straightforwardly as one tracks the original observables in terms of fermionic coherent states to the inde- 
pendent, anomalous-related, locally Euclidean coset fields a('5';'i') (I5.30ati5.30ep by differentiation with respect 



6.1 Observable quantities by differentiating with the generating source term 
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to 3l%,.^j7'^^yJ;7''^,s). This leads to the foUowing three terms (j = 1, . . . , 2iV) 



dZ[3] 



^ ^ - "-1,2 ' -^^,2 J 



a,b=l,2 



-1 



l;a6 



eq 



-27V; 



(6.1a) 

eq. (l5.22» 

(6.1b) 



-2A4.( ( i-i 7? t(S) :k-i p-i ( i + z - X ) p t-i(5) r? i-i ) (t^^i'; r, J;) 

^l,s\u,s' ' I eq. (l5.22D 



-1:12 



eq. (11211) 



-27V; 



(6.1c) 



-1:21 



2A4:( (rST(a)tf{-ip-i(i + £-if£^C-i) ^P t-i(a) r/r^ ) (Tj;'i';T^:,^) 



eq 



(6.1d) 



+ S 



Applying relations of section U] and O we can reduce and simplify the inverse of combined anomalous- and 
density-related parts in (|6.1ati6.1c|) so that one is reminded of the astonishing result that the exciton quasi- 
particles in the off-diagonal blocks of 2,''"('5'; '1') are multiplied by the one-particle operator (s)] in order 
to achieve the true exciton field 



r-l 



-1 



-l;ab 



PT-\5)r,rM (T^J;T;^;) 



J ' 3' 



(6.2) 



Ir?P 



-1 



exp{[ln(-J{),. . .]_} -i)Z 



P^p-^exp{-[...,YD]-} 



-l;ab 

p^i)(75:^';t;:^) 



3r 



In consequence, the 'Nambu'-doubled, locally Euclidean fields Aa('5';'i') have to be regarded as exciton fields 
relatively to the density-related parts. Further restriction to solely diagonal elements ais^ (t] V ; V ) of a('5'; '1') 
()5.24a|5.24bl I5.30all5r30ep allows to derive and investigate classical soliton equations from first and higher order 
variations which are determined by the prevailing coset integration measure (cf. section ri.2p . Therefore, we 
have pointed out and intensively examined the symmetries of the corresponding self-energies throughout the 
paper under inclusion of the precise time contour step ordering of the original time development operator. The 
derived equations of section 15.31 for the observables (j6.1all6.2|) can be used to study self-induced transparency 
of ultrashort coherent transients and the area theorem in combination with holography and selective Fourier 
optics for possible means of reduced absorption within matter [26]-j3Qj. Possible extensions of the presented 
approaches of sections I3.m3.3l and HHS] may include the quantization of the electromagnetic field which can also 
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6 SUMMARY AND CONCLUSION 



be represented in terms of coherent states so that corresponding HSTs lead to analogous self-energies as for 
the quartic interaction of fermionic fields. 
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A Gauge fixing for 'Nambu'-fields within the coset integration measure 

In section \^72\ with eqs. (|5.15al5.15bp we have to require that the quaternionic diagonal, anti-symmetric matrix 
elements have to vanish in the gauge combination ( (dP) P~^);i5,s;/^i,s('J']^|'; "J^'^l') ()5.15bp of the block diagonal 
matrices P^'!^^s5;iii,sii''^j^%r'^'-^ji^hi) (|4-7ll4.9fp with their contour time derivatives 'd' 

0=((ciP)p-i)" . (A.l) 

This can be accomplished by a gauge transformation (|A.2P of P('5';'i') with a quaternion diagonal matrix 
PD;/.5,s;Mi,s('^ri';'^jt£^) ()A-3IA.4p which has only non-vanishing matrix elements 9D■,^,5,s■,^Musi'^f,s '^"^l^s) / 
along the 2x2 diagonals, just in opposite to P""('5';'i') ()4.9all4.9T]) . These hermitian 2x2 matrix elements 
9D;fj.5,s;fii,s{'^l'^,s 'I'^j^s) f|A.5p aloug the main diagonal have to depend on the off-diagonal parameters of the 
ladder operators in P^^('5'; 'i'), P^'^('5'; 'i') and have to be chosen with suitable dependence in such a manner that 
the block diagonal, gauge transformed matrices iP""('5'; 'i') = P^('5';'4') P'^"('4'; 'i') (jA.2p fulfill the property 
of ( (dJ') J'^^)^5,s;^i,s('J']|s'; = (]A.ip . One has to take into account the quaternion algebra in order to 
achieve (IA.7P for diagonal elements referring to the quaternion matrix elements with anti-symmetric element 
('^2)^45/^1 of exciton quasi-particles 

paa (q-^^JsK ij^^n ') \ . ayaa (rr^'^js 'l . q-^^Ji^ \ . ( \ O') 

'^fi5,S5;f^l,Sl\''j5,Xs'-'jl,XiJ ~^ '^fJ.5,S5;fJ.l,Sl\-'j5,X5'-'ji,Xi) ' \-^-^) 

P]^,,,,s^,,,s{7'^i';7^i')=exph9j,,^^^^^^^ ; (A.3) 

9D;,„s,,us(7t:^;Tt^) = Sd(Sd('3V2');'3' / '2^(7^i';7'^i') ; (A.5) 
dT") {p-^-^""^ ('5'; '1') = '4') dP'^'i'A'; '3') F-^-^'^'i's'; '2') P^^'""('2'; '1')+ (A.6) 

+ dP?,"('5';'4') P^^''^°('4';'i') ; 

Ep-l;aa ('rr(''j5 > . a-(''i5 ' ^ J paa /'rrt'^Js ^ a-*''J5 ' ^ _ ( A 7\ 

D■,^l5,s■,^J.4,s\•'j5,X5' ■'js^xrj "■^D;fi.4,s;fii,s\-'j^,x5' j^^x^' ^-^'^ > 

/i4='eVh' 

2A'" At4,S4 

EV^/ jpaa |'n'(''i5>.T-(''i4'\ p-l;aa (cj(-^3aK cj^'^h'^ \ 

J4 = l X4 



{1 • ) 

Note that there is no summation over the contour time and coordinate space labels 'T;/^!^' and the spin label 'ss' 

in relations (IA.2|A.7p . due to the 'diagonal' property of P/>^g ^.^^ ^ ('J'^/lg ; T^^^lg ) with respect to these labels, 
apart from the electron- hole labels '^5,^*1 = 'c','h'' which restrict the degrees of freedom to the anti-symmetric, 
quaternion element {f2)^^^i- A corresponding description holds for relation (IA.6P with the quaternion diagonal 
matrix elements P''d;u,s;,i,s('^3,3 ^'^i.s )' having only the quaternion elements {f2)^,^f,^, ('f2)^2Mi' ("^2)^5/^4' ('^2)^4^1 
for respective gauge matrices P^"('5'; '4'), P^^''"'('2'; '1'), dP^^^s'; '4'), P^^''"'('4'; '1'). 
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